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Abstract 

In this paper we address the problem of representing solutions of a system of scalar 
linear partial difference equations akin to state space equations of 1-D systems theory. We 
first obtain a representation formula for a special class of autonomous systems. Then we 
show every autonomous system can be converted into the special ones by a coordinate 
transformation on Z". Using this conversion we provide representation formula for general 
autonomous systems. The representation formula we present can be viewed as multidimen¬ 
sional flow operators acting on initial conditions. These initial conditions are required to 
satisfy certain compatibility conditions. We give a full description of the set of allowable 
initial conditions. 

In our search for a general representation formula, one algebraic result plays a very 
crucial role. In this result we show that every quotient ring of the n-variable Laurent 
polynomial ring can be made a finitely generated faithful module over another Laurent 
polynomial ring of smaller dimension by doing a suitable change of coordinates. We call 
this result a discrete version of Noether’s Normalization Lemma. 


1 Introduction 

One of the most crucial discoveries in the theory of dynamical systems is the notion of state 
variables. The state variables serve, simultaneously, two purposes in the analysis and control 
of dynamical systems. On the one hand, they relate to the system’s memory, or, equivalently, 
energy storing capability, while on the other hand, they provide a recursive formula to compute 
the trajectories in the systems. The hallmark charateristic of state variables that enables the 
above-mentioned two benefits is that the representation of a dynamical system in state variables 
is a first order ordinary differential equation (ODE) of the following form: 

^ = /(x,u), y = h{x,u). (1) 

For discrete dynamical systems, that is, systems described by ordinary difference equations 
(odes), the state equation takes the form: 

x{k + 1) = f{x{k),u{k)), y{k) = h{x{k),u{k)). (2) 

Given the benefits of state variables it is desirable to have systems decribed in the form of 
equation Q. Unfortunately, often the hrst principles modelling of dynamical systems end 
up giving higher order equations |RW97] . However, for a system of ODEs, a higher order 
mathematical model can always be reduced to a state variable model (see |RW97j where this 
was proven for linear time-invariant systems). And, therefore, it is a standard practice to 
consider dynamical systems already given by a state variable model. 
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The situation drastically changes when it comes to systems described by partial differen¬ 
tial/difference equations (PDEs/pdes). For such n-D system^ as they are often called, it is 
not clear whether a state variable representation similar to equations 0 or © always exists. 
The main hindrance in obtaining such a first order representation for n-D systems is that, un¬ 
like 1-D systems, there is no fixed direction of evolution for n-D systems. This problem can be 
somewhat overcome by letting one of the independent variables play the role of ‘time’ (see, for 
example, [CZ951 [STWn2[ iNRRllj ). This approach naturally leads to an infinite dimensional 
state space |CZ95] . A serious drawback of this approach is that there are many n-D systems 
for which such a state space representation with a special variable treated as time may not 
even be possible |NRR11] . 

In this paper, we show that for every discrete n-D autonomous system, a first order repre¬ 
sentation akin to equation Q can be obtained. Our inspiration comes from the special class 
of autonomous n-D systems ~ the strongly autonomous systems. It is known that a strongly 
autonomous linear n-D system admits a first order representation of the following type: 


x{ni + 1 , ^ 2 , ...,Vn) 
x{ui,V2 + l,...,nn) 

x{ni,U2, ...,Un + l) 

w{ui,l'2, ...^Vn) 


= Aix(ni,n2,...,n„) 

= A2x(ni,n2,...,n„) 

= Aix(ni,n2,...,n„) 

= C'x(ni,n 2 ,... ,n„). 


( 3 ) 


where Ai, A 2 ,..., An are real square invertible matrices, C is a real matrix of appropriate size 
|EW89[ rFR,Z93] . Note that here, each state matrix Aj encodes the information about the evo¬ 
lution of the state variable x in the direction of the integer grid Z", over which the system’s 
trajectories evolve. The state variable x evaluated at any integer n-tuple (ni, ^ 2 ,..., n„) G Z” 
is a finite dimensional vector. So the state space here is a finite dimensional vector space. 
We show in this paper, that for a class of systems (that strictly contains the class of strongly 
autonomous systems), such a first order representation exists. However, the representation has 
a striking difference from equation (§: here the state space is usually infinite dimensional. We 
show that for these systems, there exists a non-negative integer d < n such that the state space 
is a subspace of the space of all vector valued discrete trajectories evolving over Z'’* (we call 
such trajectories d-D trajectories). And, the first order representation provides an (n — d)-D 
evolution on the state space. Interestingly, as we show in this paper, the number d is equal to 
a geometric invariant of the autonomous system: it equals the dimension of the characteristic 
variety of the system. As a result of this, the representation of equation ^ turns out to 
be a special case, because the defining property of strongly autonomous systems is that the 
dimension of their characteristic varieties is zero, therefore, d = 0 in this case, and so the state 
space turns out to be finite dimensional. 

We call this special class of autonomous systems strongly relevant of order d. The char¬ 
acterizing property of this class is a certain algebraic property. The crucial benefit of this 
algebraic characterization is that every discrete n-D autonomous system can be transformed 
in to an ‘equivalent’ strongly relevant system by a change of coordinates in Z*^. This enables 
us to provide a first order representation for any discrete autonomous n-D system. A crucial 
role in this transformation is played by a normalization process, which we call the discrete 
Noether’s normalization lemma (DNNL), for it is an analogue, applied to Laurent polynomial 
rings, of the classical Noether’s normalization lemma, which is applied to polynomial rings. 

The above-mentioned first order representations can be utilized to provide explicit solution 
formulae for a given discrete autonomous n-D system. We provide these formulae as algorithms 
for solving system of autonomous pdes in this paper. Interestingly, a parallel can be drawn 
between these representation formulae and the classical result of Ehrenpreis-Palamodov integral 

^This way, systems described by ODEs/odes are called 1-D systems. 
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representation theorem for continuous autonomous n-D systems |Bjo79 , lObeQOi [Stu02| . The 
integral representation formula can be viewed as an evolutionary equation of initial conditions 
specihed on the characteristic variety (which is a d-D object in C”). The representation 
formulae of this paper, too, are evolutionary equations on state spaces, which are subspaces of 
d-D trajectories. However, a striking difference is that, in our case, since the state space is just 
a subspace of d-D trajectories, the initial conditions are specified over as opposed to the 
ones specified over the characteristic variety in the integral representation theorem. Further, 
our evolutionary equation does not involve an explicit integration. These two facts make our 
representation formulae easily implentable on a computer. 

The resuts presented in this paper are generalizations of those in [PP13j . which deals with 
the special case of n = 2. It is well-known that often such an extension from n = 2 to a 
general n is not straight-forward, owing to the fact that dimension 2 algebras have certain 
special features that do not carry over to algebras with higher dimensions (see |PS98j for some 
of these special features of 2-variable polynomial rings). In our case, the main difficulty arises 
in obtaining the DNNL for higher dimensional rings. Another crucial difference between 2-D 
and n-D is that, in case of general n, the number d can be anything between 0 and n — 1, while 
in case of 2-D the number d can only be 0 or 1. 

The paper is organized in the following manner: in the next section (Section]^ we discuss 
the notation used in this paper, and provide some preliminaries required for the main subject 
matter of the paper. Then, in Section we consider a special class of discrete autonomous 
n-D systems and show how a first order representation of such systems can be obtained. Using 
this representation we then provide an explicit solution formula, alongwith an algorithm, to 
compute the trajectories in such systems. In Sectionj^we state and prove the discrete Noether’s 
normalization lemma. Then using this normalization we provide, in Section a first order 
representation for general autonomous discrete n-D systems. As a consequence, we provide an 
explicit solution formula for trajectories in general autonomous discrete n-D systems, and also 
an algorithm to compute them. Finally, in Section]^ we provide some interesting insights on 
the set of allowable initial conditions that are used in the solution formulae. 


2 Notation and preliminaries 

2.1 Notation 

We use M and C to denote the helds of real and complex numbers, respectively. Consequently, 
C” denote the n-dimensional vector spaces over M and C, respectively. The set of integers 
is denoted by Z, and denotes the set of n tuples of elements in Z. In this paper, our main 
object of study is a particular class of multiply indexed sequences of elements in for some 
positive integer q. We denote the set of sequences in M indexed by Z” by the symbol W, i.e., 
W := {tc : Z” —)■ M}. To denote the set of vector valued sequences indexed by Z"', we use 
i.e., W'^ := {w : Z"" —)• M'^}. The Laurent polynomial ring in n indeterminates cti,(T 2 , ... ,(7^, 
usually written as will be denoted, in this paper, by A, and the same 

in the first d indeterminates cri,a 2 ,. ■., cr^, written as M[(Tf ..., will be denoted by 

Ad- We use A‘^ to denote the free module of rank q over A, where the elements of A'^ are 
g-tuples written as row-vectors. For a set S, we use to denote the set of (m x n) matrices 

with entries from the set S. We often use col(») to stack up entries one above the other to 
make a column vector. The single letter cr is often used to denote the n tuple (cii, cj 2 , ..., cr„). 
Further, for an integer n tuple v = (z^i, 1 / 2 ,..., r'n) £ the symbol denotes the monomial 
this paper, we follow the bar notation to denote equivalence classes: for 
r(<T) G A^ and a submodule TZ C A^^, we use r(cr) to denote the equivalence class of r{a) in 
the quotient module A^ jTZ. 
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2.2 Discrete n-D systems 


By discrete n-D systems, in this paper, we mean systems described by a set of n-D linear 
partial difference equations with constant real coefficients. Such partial difference equations 
are often described using the n-D shift operators fii, cj 2 , ..., <7^. These shift operators act on 
the n-D discrete ‘trajectories’ (i.e., real valued sequences indexed by Z”) w G as follows: 
for n' := (n^, ^ 2 ,... v := (ni, ^ 2 ,..., n„) G Z"- 

vS){y) = w{y ^ v'), (4) 

where the sum v^v' is the standard element-wise sum: v = (ni -|- ^2 -|- ^ 2 ,..., -|- n^). 

This definition can be extended naturally to define the action of A, the Laurent polynomial 
ring in the shifts, on W. And, likewise, a row-vector r{cr) = [ri(cr), r 2 (cr), • • • ,rq{a)] and a 
column-vector w = col(rci, t(; 2 ,... ,Wq) G we define r{a)w := This defines 

the action of the row module A^ on 

In this paper, as a convention, we always consider a vector valued trajectory as a column 
vector. Consequently, any tuple of difference operators that act on a column vector of trajecto¬ 
ries and produces scalar trajectories, is written as a row, so that this action can be compactly 
written as a ‘dot’ product of a tuple of difference operators and that of trajectories. So, as a 
convention, in this paper, elements from a module over an operator ring A (or Ad, as the case 
might be) are assumed to be row vectors. 


2.3 The kernel representation 

Following Willems’ 1-D terminology [Wil91| . we call the collection of trajectories w G that 
satisfy a given set of partial difference equations the behavior of the system, and denote it 
by *B. The above description of the action of A^ on gives the following representation of 
behaviors of discrete n-D systems: 

^ ■= {w \ R{<t)w = 0}, (5) 

where R{(t) G A*^’^. The above equation ([^ is called a kernel representation of IB and written 
as IB = ker(i?(cr)). Note that various different matrices can have the same kernel. Importantly, 
all matrices having the same row-span over A result in the same behavior. This leads to the 
following equivalent definition of behaviors: let R{cr) G A*^^ and TZ := rowspan (R{cr)), 

iB(7^) := {w G I r{a)w = 0 for all r(a) G TZ}. (6) 


The submodule TZ generated by the rows of a kernel representation matrix is called the equation 
module of *B. It was shown in |Obe90| that the submodules of A'^ and discrete n-D behaviors 
that have q number of dependent variables (that is w{i') G for all n G Z") are in an inclusion 
reversing one-to-one correspondence with each other. We denote the set of all discrete n-D 
behaviors with q number of dependent variables by 

In mathematical modelling of systems, it often turns out to be beneficial to incorporate 
auxiliary dependent variables called latent variables |Wi191) . A representation of the behavior 
involving both the variables of interest (the manifest variables |Wil91| i and the latent variables 
is called a latent variable representation. The following equation is one such representation: 


IB 


= |t(; G 3 £ G such that R{a)w + M{a)l = o| 


Note that the state variables are a special type of latent variables. 


(7) 
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2.4 Autonomous systems 

In this paper, we deal with a special type of discrete n-D systems, namely autonomous systems. 
Among several equivalent definitions of n-D autonomous systems (see |FI1,Z931 IPS981 rVaini| h 
in this paper, we stick to the algebraic definition of |PS98j . which we state below as Definition 
In Definition and later we need the notion of characteristic ideal of a discrete n-D system. 
Let the behavior of a discrete n-D system be given by a kernel representation OS = ker(i2(cr)) 
with R{cr) G The characteristic ideal of 55, denoted by X(iB), is defined as the ideal of 

A generated by the {q x q) minors of R{a). If the number of rows in R{o’) is less than the 
number of columns g, then X(iB) is defined to be the zero ideal. 

Definition 1. A discrete n-D system is said to be autonomous if the characteristic ideal I{^) 
is nonzero. Further, an autonomous behavior is said to be strongly autonomous if the quotient 
ring A/X(*B) is a finite dimensional vector space over M. 

2.5 The quotient module 

Given a behavior 55 = ker(i?(cr)), let TZ be the submodule of A'^ spanned by the rows of R{cr). 
We define 

M :=AyTZ, 

and call it the quotient module of *B. This quotient module AA plays a central role in this 
paper. We often let elements from AA act on iB. This action is defined as follows: for m € AA, 
the action of m on re G 55 is defined to be the action of a lift of m in A'^ on w. For example, 
let ^(cr) G A'^ be such that r{cr) = m G AA, then 

mw := r{a)w. (8) 

Note that m may have several distinct lifts in A^, but all of them have the same action on 
tc G 55. This is because two lifts differ by an element in the equation module TZ. Therefore, 
the action of this difference on any tc G 55 produces the zero trajectory, for the action of any 
element in TZ is zero on 55. Thus, the definition of action of AI on 55 is well-defined. 

Now note that it follows from Definition above that is autonomous if and only if the 
quotient module Al is a torsion module, i.e., for every m[a-) G AA there exists a /(cr) G A 
such that f{cr)m{cr) = 0 G Al. In that case we get the following ideal called the annihilator 
ideal of AA. 

ann (Al) := {/(cr) G A \ f{cr)m = 0 for all m G Al}. 

3 A special case: A/ is a finitely generated module over a 
smaller ring 

We first consider a special case. Let the n-tuple of shift operators be partitioned into two 
tuples: cr = (cri,cr2), where <ti = (cti, (T2, ..., Urf) and <T2 = {ad+i,crd+2, ■ ■ ■ ,crn)- In this 
section, we consider the case when the quotient module Al happens to be finitely generated as 
a module over the subring M[cri,<Ti“^] (that is, the Laurent polynomial ring generated by the 
d-tuple of shift operators cri). We denote this smaller ring by Ad. 

The motivation for considering this special case comes from the intention to mimic the 
construction of states for strongly autonomous systems. Recall that, for strongly autonomous 
systems, the quotient module Al happens to be a finite dimensional vector space over M. It is 
already well-known that, in that case, the n-D behavior admits a state-space representation of 
the form given by equation ([^ (see []). However, there are many n-D autonomous systems that 
are not strongly autonomous. For such systems the quotient module is an infinite dimensional 
vector space over M. Among these systems, there is a special class, for which Al is a finitely 
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generated module over a smaller ring Ad- We formally define such systems as strongly relevant 
system of order d. 

Definition 2. Let 55 G he an autonomous system. Then 55 is said to be strongly relevant 
of order d if the the quotient module A4 of ^ is a finitely generated module over Ad- 


For example, the following scalar 2-D system 


55 = ker 


(t| - 2(T2 + 1 
(T 1 CJ 2 — (T2 — CJi + 1 


is a strongly relevant system of order 1. On the other hand, the scalar 2-D system 


55 = ker [(Ticj 2 — (T 2 — cji -|- 1] 

is not strongly relevant of any order. Note that a strongly relevant system of order 0 is nothing 
but a strongly autonomous system. 

The assumption that is a finitely generated module over Ad gurantees the existence of 
certain matrices which will be crucially used in obtaining the desired representation formulae. 
We describe the construction of these matrices and their properties below. 


3.1 The matrix of relations, X(cri) 


The quotient module At being finitely generated over Ad implies that we can find a finite 
collection of elements from At, say {gi,... ,gr}., that generate At as an A^-module. This 
choice of generators gives us a surjective A(i-module morphism 


V’ : ^ At 

ei 1 -^ Qi 


(9) 


for 1 ^ i ^ r, where is the standard basis row-vector 

Ci = [0 0 • • • 1 • • • 0 ] G AJ;. 

t 

i^^ position 

Clearly, ker('0) is a submodule of AJj. Since A^ is a Noetherian module, this submodule is 
finitely generated. Therefore, we can find a matrix X{cr-j_) G A^^^ whose rows generate ker('0) 
as an A^-module. Thus At turns out to be isomorphic as an A^-module with the quotient 
module A^/rowspan(A(<Ti)): 

AJj/rowspan(X(<Ti)) ~ At as A(i-modules. (10) 


We call A(<ti) a matrix of relations of the generators of At as an A^-module. 


3.2 Matrices representing the (n — (i)-tnple of shift operators ad+i,... ,an 
Now consider the A^-module morphisms defined for 1 ^ j ^ (n — d) as 

gj : M —^ At 

m I—)■ fj,j (m) := (Jd+j • (11) 

These maps are clearly A^-module morphisms. Since At has been assumed to be finitely 
generated over Ad, in order to specify these maps it is enough to specify their actions on the 
generators {gi,... ,gr}. Suppose that for 1 ^ f ^ r and 1 ^ ^ (n — d) 

Tj (di) — ®il,j(^l)dl T ®i2,j(^l)d2 T ■ ■ ■ T 0,ir,ji^l)9r, (t^) 


6 



where ajij(cri), Oi 2 ,j(cri),.. 

• ; O^irJ 1) 

G Ad- We denote 

by Aj{ai) 

the following matrix 


ttiij(cri) 

oi2j(o'i) ••• 



Aj{cTi) = 


a22,jr(o'i) ••• 


G A'd^F 



®r2,j (^l) 

CLrrji^l^ 



(13) 


Equation (12) can now be written collectively for all i as 



f 

’ 91 ' 

\ 


’ 91 ’ 

9j 


92 


= ^i(o-l) 

92 


V 

. 9r _ 

) 


. 9r _ 


(14) 


It then follows that multiplication by cr^+j is represented by the matrix Aj {a\). We prove this 
in the following lemma. 

Lemma 3. Suppose A4 is a finitely generated module over Ad- Let {gi,... ,gr} ’L Ai be a 
generating set for Ai as an Ad-module. Consider : A^ AA as given in equation p). 
Further, for 1 ^ j ^ [n — d), let p,j and Aj(cri) be as defined by equations 0 and 
respectively. Then the following diagram commutes. 

V’ 

^d - M 

Aj{cri) 4, I pLj (15) 

Ip 

^d - M 

where Aj{a\) : A^^ —)• A^^ is defined by right-multiplication by the matrix Aj{a\) to the row- 
vectors in A^d- 

Proof : For any r{a\) = [ ri(<Ti) r 2 ((Ti) • • • rr(<Ti) ] G Ad we want to show 


From the definition of the map it follows that 


'(/;(r(cri)) = rfiafijgi + r2{cTfi)g2 H-h rr{(Tfi)gr. 


(16) 


Operating pLj on both sides of equation (16) and then utilizing equation (14) we get 
/ij(V^(r(cri))) = /ij (ri(<Ti)fi(i+r2(<Ti)fi(2 H- ^ rfia^^gr) 




f 

’ 91 ' 

\ 

[ ri(cri) r2(<Ti) •• 

• Vrierfi) ] pij 


92 




V 

. 9r _ 

/ 


= [ n(cri) r2(cri) ••• r^(cri) ] ^^(cri) 


91 

92 

9r 


= 'fiir{o-i)Aj{ai)). 


□ 


The matrices Ai[ax), A 2 {(Tx), ..., ^n-d('^i) called companion matrices. In the lemma 
below we list out some useful properties of these companion matrices. 
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Lemma 4. Consider if; and Hj, Aj{ai), for 1 ^ j ^ {n — d), as defined by equations 
and (13), respectively. Then we have the following: 


1. The matrices Ai{a\), A 2 {a--i),... ,An-di(^i) commute pairwise. 

2. There exists a set of generators {gi,g 2 , ■ ■ ■ ,gr} for A4 as a module over Ad such that 

each of the corresponding Aj((Ti) ’s is invertible in (that is, unimodular over Ad). 

3. With each of Aj{a-i)’s invertible, we have the following commutative diagram for any 

e 




•^d 


a: 


M 
i « 
M 


Uo 

°h2 ° 


yn — d\ 

° l^n-d ) 


(17) 


4 . With each of Aj{a\) ’s invertible, for every {vi, U 2 ,..., Vn-d) £ , there exists a matrix 

E{a-i) G A^^^ such that 


X{<Tl) 



EMXM. 


(18) 


Due to its length and technical nature we have given the proof of this lemma in the ap¬ 
pendix. 


3.3 The output matrix, C'(cri) 


The output matrix is related with the image of the q x q identity matrix under the canonical 
map A^ ^ M.. We elaborate on this now. Let s* denote the ith standard basis (row-)vector 
for A'^, and si its image under the canonical map A'^ Ai. Now, consider again the map 
if : A(i -P- M. given by equation g. Since if is surjective, for each Si G A^, with 1 ^ i ^ q, 
there exists Ci(cri) G A(i such that 

fi’{ci{o-i)) = Xi. 


The output matrix C{a\) G is obtained by stacking these Ci(cri)’s one below the other. 

That is. 


C{CT1) 


ci(cri)' 

C2(o'i) 




(19) 


OqM. 


3.4 The first order equation 


The key observation to obtaining a representation formula for strongly relevant systems (that 
is, M. being finitely generated over Ad) is the fact that such systems admit a first order 
representation. This fact is very much akin to the well-known notion of state space in ID 
systems theory. However, unlike ID systems, the state space in this case is often not a finite 
dimensional vector space. In this subsection, we show the construction of the state space and 
the corresponding first order representation of the original system, using the three type of 


matrices introduced in Subsections 3.1, 3.2 and 3.3 










Recall that for strongly relevant systems, the quotient module is a finitely gener¬ 
ated A^-module. This means we have a matrix of relations the companion matri¬ 
ces Ai(cri), A 2 (cri),..., and the output matrix C'(cri). Define the following d-D 


behavior 


3^ := <^ X G 


X{cri)x = o| 


We call ^ the state space. 

Now note that, as a consequence of statement 4 of Lemma[^ we have, for all 1 ^ j ^ {n — d) 


Ai{a\)^ C ^. 


( 20 ) 


That is, is invariant under the action of each of the companion matrices ^ n-d- 

Therefore, each Ai[a\) can be thought as an endomorphism of ^. We are going to use these 
companion matrices to dehne a first order (n — (i)-D system over For this purpose, let us 
denote by x an (n — (i)-D trajectory over Alternatively, x is a map from —>■ S'. We 
now describe the first order system over S' as given below: for (z^i, z^ 2 , • ■ ■, ^n-d) £ Z”'"'^, 


x(l^l 1 , 1 / 2 , 

■ ■ ■ 5 d) 

= Ai(<Ti)x(z/i, 1 / 2 ,.. 

• 1 ^n—d) 

x(l/i,Z/2 -h 1, 

• • ■ ; ^n—d) 

= A2(cri)x(z/1,1/2,.. 

• ? ^n—d) 

x(z/i,Z/2, . . . , 

^n—d “1“ 1) 

= A„_rf(cri)x(j/i,z/ 2 ,. 

• • ? ^n—d 


Remark 5. Note that x can also be identified as an n-D trajectory. This identification is 
done as follows. Suppose v = {ui, 1 ^ 2 , ■ ■ ■, ^'n) £ is arbitrary. Then x, identified as an n-D 
trajectory, maps this u to {x{i'd+i, i^d+ 2 , ■ ■ ■, 1 ^ 2 , ■ ■ ■, ^d)- In what follows we do this 

identification implicitly throughout the rest of the paper. 

With the identification mentioned in Remark]^ it follows that (^(crijx can be viewed as 
an element of We summarize these constructions in Table 1 below. 


TABLE 1 


s: := ■ 

X G 



A(cri)x = 0 



X 

z"-' 

^ ^ s:. 


For (i/i, z/2,..., Vn-d) G 

x(l^l -M,i/ 2, ... , r'n-d) = Ai(cri)x(z^i,l/2,---,t'n-d) 

x(z/i,z /2 + 1 ,.. . , r'n-d) = A2(<Ti)x(z/i,ZZ2,...,r'„_rf) 


x(i/i,Z/2, . . . ,^'n-d + 1 ) = A„_rf((Ti)x(;zi,zy2, . . . ,Z/n-d)- 

w = C(cri)x. 

With this first order system of equations we are now in a position to state and prove the 
representation theorem for the strongly relevant systems of order d. 


Theorem 6. Suppose ^ ^ 23 is a strongly relevant autonomous system of order d. That is, 
the quotient module M. of is finitely generated as a module over Ad. Then there exist a 
positive integer r and the following matrices: 

. A(^i) G 


• Ai(cri), ^2(0-1), . . .,An-di(Tl) G A"'/'', 

. C(^i) G Af ^ 

such that IB admits the following representation: 


= 


{rcG 


3 xG 


such that equations in Table 1 are satisfied^ (21) 
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We need the following lemma to prove Theorem The lemma concerns a test, involving 
the matrices X(cri), ^i(cri),..., C'(cri), as to whether a given vector r(cr) G A'^ 

maps to 0 under the canonical surjection ^ Ai. 


Lemma 7. Let Ai = A'^jTZ be a finitely generated module over Ad with {gi,... as a 
generating set. Let : A^ Ai be given by equation ^ and let X{a-i) G A^^'^ he a matrix 
of relations of {gi,... ,gr}- Further, for 1 ^ j ^ (n — d), let Aj{(T\) he as defined by equation 
(13), and let C{(J\) G A^^'^ he as defined in equation (19). Suppose f{(j) G A^ is given by 


/(cr) = 




a,,cr 


with au G Mrow and only finitely many of {a^} nonzero. Then the following statements are 
equivalent: 

1 . f{a)en. 

nil ITjZt G ker(V^). 

3. There exists a matrix F{a\) G A^^^^ such that 


d n—d 

Aj{ary<^+^ = F(cri)X(cri). 

i=l 1 


Proof : The proof follows by noting how a typical element, say /(<t) g A'^, corresponds to 
a vector in A(i, such that the images of these two elements under the maps A'^ Ai and 
Ai, respectively, coincide with each other. That is, given /(<t) g A^, how we can 
find out an r(<Ti) G A(i such that 

^(r(cri)) = f{a). 


To get this correspondence, first consider /(cr) to be a monomial, that is, /(cr) = aa’-', where 
a G Mfow and v G TP. In this case, the image of /(cr) under the canonical map A^ -» is 


given by 

/(cr) = a^. (22) 


Let us denote the component v by Vi. 
rewrite equation (22) as 


Then by partitioning cr into cri and (T 2 we can 


/(cr) = aa^ 

d n 

= n 

i=l i=d-\-l 

n 

= n (23) 

i=d-\-l 


where ui denotes {ui, U 2 ,..., Vd)- Now note that for a G M?ow, it follows form the definition of 
the output matrix C(cri) that 

a = fi{aC{a--i)). 

Therefore, since ip is an A^-module morphism, we also get 

acTi^^^ = p {ao'x''^C{ai)). 
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Putting this in equation (23) we get 


/(cr) = cFi’^^ 

n 

i=d-\-l 

= o o • • • o ^^n-d) ^ (ao-i‘"iC'(cri)) 

/ n—d \ 

= V' facTi'^i(7(0-1) , 

where the last equality follows from statement 3 of Lemma 

For the case when f{cr) G A'^ has more than one terms, we note that ijj is Tl^-linear. 
Therefore, the above argument can be applied term-wise. Thus, for an f[cr) G A^^ given by 

/(cr) = X] 


we get 


d n—d 

/(cr) = V' I X] ^^W^TC{(Ti) Aj{aiY'i+i 

i=l j=l 


(24) 


The result of the lemma then easily follows by noting that f{cr) G 7^ if and only if /(o-) = 
0 G At. □ 


With Lemma we are now in a position to prove Theorem 
Proof of Theorem^^- Suppose 53 G is an autonomous behavior with quotient module At. 
Let us define 5Saux as the set of trajectories admitting the representation given by Table 1. 
That is. 


:= 


|t(; G ^ 3 X G ^ such that equations in Table 1 are satisfied| . (25) 


Then we have to show the following equality of two sets of trajectories: 


= 53 


aux* 


(To show C 53aux) We have to show that tc G 53 implies that there exists x G (M^")^ 
such that equations in Table 1 are satisfied. For the sake of easy referencing we reproduce the 
equations given in Table 1: for all {ui, 1 / 2 , ■ ■ ■, ^n-d) £ x satisfies 



x(z^l,Z/2,. 

..., Vn-d) G and 

(26) 

x(l^l 1, 1/2, 

■ ■ • ? ^n—d) 

= Ai(<ti)x(i/i, 1 / 2 , ■ ■ 

• ? ^n—d) 1 

x(z/i,Z/2 + 1, 

• ■ • ? ^n—d) 

= A2(<Ti)x(i/i, J/2, .. 

(27) 

x(z/i,Z/2, . . . , 

^n—d “1“ 1) 

= A„_rf(cri)x(j/i,z/ 2 ,. 

• • , ^n— d), J 


with 

w = (7(cri)x. (28) 

In order to get this x suppose { 51 , 52 ) • • • > 5r} C Af is a generating set for Af as a module over 
Ad- Now, let us define for re G 53 the following r-tuple of n-D trajectories: 


51 

52 


_ 9r _ 


(29) 
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We claim that the x defined above satisfies equations (26), (^27^ and 


We first show that x(z^i, U 2 , ■ ■ ■, Vn-d) £ ^ for all (z^i, 1 / 2 ,..., Vn-d) £ Note that it 


follows from the defining equation (29) of x that 


x(i^i,z^2, • • ■,J^n-d) = 


91 

92 

9r 


w{»,Ui,U2,...,12n-d)- 


We make the matrix of relations X(cri) act on x(i/i, 122 , ■ ■ ■, I'n-d) to get 


X{a-i)y.{vi,V2,...,Vn-d) = X{ai) 


= X{a^) 


91 

92 

. 9r _ 

91 

92 

9r 


W{*, 121,122, . . ■,12n-d) 


W{», 121,122, . . .,I2n-d) 


= 'ijj {X{cr-l_))w {•, 121 , 122 , ... , 12 n-d) 

= 0 , 

because rowspan(X(cri)) = ker{ijj). Therefore, x satisfies equation 


We now show that x satisfies the system of equations (27). Let i G {l,2,...,(n — d)} and 
{ 121 , 122 ,... ,i2n-d) £ be arbitrary. Recall how x is treated as an element in see 

Remark Under this identification, it follows that 


x{l2i, 122,. . . ,12i + l,. . . , Vn-d) = CFd+iX{vi, 122,..., 12n-d) ■ 


(30) 


Applying the defining equation (29) of x to equation rt30| above we get 


X{l2i,122,...,12i + l,...,12n-d) = 


ad+iX{l 2 i, 122 , . . .,12n-d) 

91 






( 


— 9‘i 


92 

. 9r _ 

91 

92 

9r 

91 

92 


w{», 121,122,..., Vn-d) 


w{»,Vi,V2, . . . ,Vn-d) 


\ 


V L . 

= Ai(cri) 


w{»,Vi,V 2 ,...,Vn-d) 


91 

92 


9r 


w{»,Vi,V2, . . .,Vn-d), 
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where the last equation follows from equation (14). Therefore, we get 

x(zyi, 1/2, . . . , I/j + 1, . . . , I'n-d) = Ai(<Ti)x(l/i, Z/ 2 , . . . , Z/j, . . . , Vn-d) 

for all 1 ^ ^ (n — s) and [vi, z/ 2 ,..., Vn-d) £ That is, x satisfies equation (27). 

Finally, we have to show that x and w satisfy equation (28). This follows from the definition 
of the output matrix C'(cri). Multiplying both sides of equation (29) by C'(cri) we get 

C'(cri)x = C(cri) w 


9r 

= 'll) {C{ai)) W = IqW = w. 

This concludes the 53 C 5Saux- 

(To show 53 ^ ®aux) We now show that if we obtain a rc by solving the system of 
equations in Table 1, that is, equations (26), (27), (28), then that rc is in !B. Let 77 C 
be the equation module of 53. Then, it is enough to show that the action of every element 
r(cr) G 77 on any w G 5Baux produces the zero trajectory. 

Now, let r(cr) be any element in 77 given in an expanded form as 


r(cr) = ^ aucr 
i/eZ" 


1/ 


where au G Krow and all of them but finitely many are zero. Once again, we partition the shift 
operators into two, namely {ai ,..., ad} and {ud+i,..., an}, to rewrite r(cr) as 


n—d 


(a) =X] n n 

i=l 2=1 


(J 


^d+i 

d+i ’ 


where (z/i, z/ 2 ,..., Un) = Then, this r(cr) is made to act upon an arbitrary w G 53aux to get 
the following: 


d n—d 


riCTjW = 


a 


d-\-i 


w 


E 

2 = 1 2 = 1 / 

^ d n—d \ 

E n ^ C{ai)^ 


a 


^d+i 

d-\-i 


X 


VZ/'GZ^ 2 = 1 2 = 1 / 

^ d n—d 

E n 

V2^GZ^ 2 = 1 2=1 / 

^ d /n—d \ 

E n 

i=\ Vi=l J 


(31) 


Before proceeding, note that 


= ^i(cri) 


'X, 


because x satisfies the first order system equations ( |27[ ). Putting this in equation (31) above 
we get 


/n—d 


r{a)w = I ^ a^]J<'C'(cri) I ^i(cri)^''+‘x 
Vz^GZ’" 2=1 


\ 2=1 
n—d 


E n X. 


(32) 


\ueZ" i=i 


i=l 
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Since r(cr) G 7^, by Lemmaj^ there exists F{cri) G such that 

d n—d 

^ C(<Ti) ]J Ai{aiY'^+- = F(cri)X(cri). 

t^GZ^ 2 = 1 2=1 


Making use of this observation in equation (32) above we get that 

( d n—d \ 

ai.]^crf C'(cri) | x 

Z^GZ'^ 2=1 2=1 / 

= F(cri)X(cri)x = 0, 


where the last equation follows from the fact that x satisfies equation (26), that is, 


x(z/i,..., Un-d) e ^ = ker(X(<Ti)) 


for all (z^i,..., Un-d) £ Therefore, rc G *B, and this concludes the 53 5 5Saux part of the 

proof. □ 


We illustrate this result in the following example. 


Example 8 . Let us consider the 3-D scalar system given by the following kernel representation: 


53 = ker 


(Tg - 2(73 -I- 1 
(t| — 2(72 + 1 

(7i(T 3 — (7i — 72 — <73 -|- 2 


In this case, clearly, A = M[ 7 i, 72 , 7 ^^, 73 , 7 ^^]. The role of the equation module TZ is 

played by the ideal a = (73 — 273 -|- 1 , 72 — 272 -|-1, 7173 — 7 i — 72 — 73 -|- 2). It can be checked that 
the quotient ring A4 = Aja is finitely generated as a module over = M[ 7 i, 7 ( 1 “^]. Therefore, 
53 is strongly relevant of order 1 . Let us compute the required matrices. As in the main text, 
let Ai(7i), ^ 2 ( 72 ) represent multiplications by 7i, 72 , respectively, in M.. First, we need to hx 
a generating set for ^/a as a module over ^ 1 . The set {1, 7 ^, 7 ^, 0 ^ 73 } meets the requirement 
for an appropriate generating set. The matrices Ai( 7 i) and ^ 2 ( 71 ), in this case, are given by 


Ai(7i) 


■ 0 

1 

0 

O' 


■ 0 

0 

1 

O' 

-1 

2 

0 

0 

, ^2(0-1) = 

0 

0 

0 

1 

0 

0 

0 

1 

-1 

0 

2 

0 

. 0 

0 

-1 

2. 


. 0 

-1 

0 

2. 


Note that both Ai( 7 i) and ^ 2 ( 71 ) are invertible over .Ai. A matrix of relations, A( 71 ), for 
the chosen generating set is given by 


A(7i) = 


I -1-11 
— 7l -|-2 —1 7l — 1 0 


It is easily checked that statements 3 and 4 of Lemma are satisfied here, that is. 


1 O' 
71-1 1 


A(7i)Ai(7i) = 

The state space 2^ is thus given by 

3F = ker A( 71 ) C 


A(7i), A(7i)Ai(7i) = 


1 0 

-1 1 


A(7i). 
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The output matrix C'(cJi) turns out to be 


C{ai) = [l 0 0 O] . 


Now, according to Theorem tc G 55 if and only if there exists x 




^ such that 


x(z^i +1,1^2) = Ai{ai)yL{vi,V 2 ) 
x(z/i,J/2 + l) = A2{(Ji)y.{vi,V2), 


(33) 


and w = C'(o'i)x. 

Remark 9. It is interesting to note that in Example by utilizing the first order represen¬ 
tation, an explicit solution formula for tc G 55 can be given. Note that every x : Z^ —>■ JE” 
satisfying equation (33) can be obtained as 


x(i^i,z^2) = Ai{aiY^A2{criY'^x, 

where x G (M^)^ satisfies X{ai)x = 0, that is, x G Now note that w = C{ai)x. Therefore, 
by the identifiaction as per Remark 


w{vi,V2,vYj = (C(cTi)x(i/2,Z^3))(i^i)- 


With this, now we can make use of the explicit solution formula for x to obtain a similar 
solution formula for w. This is given by 


w 


= [1 000] 


■ 0 

1 

0 

O' 

^2 

■ 0 

0 

1 

O' 

-1 

2 

0 

0 


0 

0 

0 

1 

0 

0 

0 

1 


-1 

0 

2 

0 

. 0 

0 

-1 

2. 


. 0 

-1 

0 

2. 




x{vi) 


(34) 


We utilize this solution formula of equation (34) in Algorithm 11 below. 


As a corollary to Theorem we get an equivalent latent-variable representation for every 
strongly relevant autonomous system of order d. 


Corollary 10. Let^ G be an autonomous system that is strongly relevant of order d. Then 
there exist 


• a positive integer r, 


• a matrix X (cri ) G A*/"', 

• a set of n — d square matrices Ai(cri),..., A„_rf(cri ) G A'X'', 

• a matrix C(<ti) g A'Y'^ 

such that 55 admits the following latent variable representation: 


/ 


r ^(‘^i) 

0 ' 


' 




0 



n;G 

3 X G satisfying 

(^d+2lr — A2{(Ji) 

0 


X 

W 

= 0 



^nlr 

0 





[ -CicTl) 



> 


(35) 


15 



















3.5 An algorithm for obtaining all solutions 

Theorem provides a useful algorithm for obtaining all solutions to a strongly relevant au¬ 
tonomous system. We describe this algorithm now. The algorithm involves two steps, which 
we call regularization and recursion. 

Given a strongly relevant autonomous system of order d, we first obtain the matrices 
A(cri), Ai{a\), ..., A„_rf(cri), and C[a--i). We call this process the regularization step. This 
step requires finding out a suitable set of generators for the finitely generated A^-module Ai. 
One has to be careful in choosing the generating set so as to ensure that the companion matrices 
Ai(cri),... ,An-d{o-i) are all invertible over Ad- This can be implemented on a computer by 
standard computer algebra packages. In this paper, we do not go into the details of issues 
related to such an implementation. 

The regularization step leads to the recursion step, where the companion matrices are used 
recursively on trajectories in = ker(A(cri)) to obtain solutions to the original system. The 
recursion indeed gives the solutions because of the first order system equations listed in Table 
1 . 

In order to see that the recursion gives an explicit solution formula for trajectories in the 
originial system, first note that solutions to the first order system defined on by equations 
given in Table 1 can be written explicitly as follows: 

n—d 

x(l/i, Z^2, • • • , l^n-d) = n (36) 

i=l 

for arbitrary ... ,n^-d) G where x G satisfies A(cri)x = 0. Since w = 

C(<Ti)x, recalling the identification of x as a trajectory in we get 


w{vi,n2, ...,Un) = {C{(Ti)x{ud+l,l'd+2, • ■ ■ , ■ ■ ■ , Vd), 


for [ui, U2,..., Vn) G Z”. Utilizing equation (36) now we get 


/ n—d \ 

w{vi, V2,...,l'n)= ( C'(cri) Ai{axY‘^+'^x I (z^i, 1/2, ■ ■ • , Vd). 


(37) 


Algorithm |11| below elaborates how this equation (37) can be used to obtain all solutions to a 
given strongly relevant autonomous system. 

Algorithm 11. (Solution formula for strongly relevant systems of order d) 
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Regularization 

Input: An autonomous system = ker(i?(<T)), where R{cr) G 

Computation: 

Obtain a generating set G := {< 71 ,..., Qr} C Ni. 

Obtain the matrix of relations X{ai) G . 

Obtain the companion matrices Afla-fl ),..., A„_d(<Ti) G A^^'' 

(ensure that these matrices are invertible in .4^). 

Obtain the output matrix C'(cri) G A^^^. 

Output: The matrices X((Ti), Aflcrfl),..., An-d{(^i)-, and C{(t-i). 

Construction of iX' 

Input: The matrix X{a\). 

Computation: Solve X{a-i)x = 0 for x G . 

Output: The set = {x G X{a-i)x = 0}. 

Recursion 

Input: An x G ^, and the matrices Aflcrfl ),..., A„d(<Ti), and C{(j{). 

Computation: 

for V = {vi, 1 ^ 2 ,..., Vn) G 1A set 

w{u) = (^ 0 ( 0 - 1 ) Ai{cr-i fl‘‘+ix'^ i’^ 1 , 1 ^ 2 , ■■■,i^d)- 

Output: The sequence w := {w{u)}y^z'^. 


Remark 12. From the algorithm above, it makes sense to call the d-D trajectories x G JE” 
as initial conditions. Indeed, the companion matrices act on x as a (n — d)-D flow operator 
to produce all the solutions in the original n-D system 55. We call this condition that x € ^ 
that is X{cri)x = 0 as compatibility condition. Note that ideally the compatibility condition 
must be satisfied at every step of the recursion (see Theorem]^. However, it turns out that 
checking for compatibility only once - at the beginning of the recursion - is enough. This is 
because the set of trajectories satisfying the compatibility conditions, that is, iX" = ker(X(<Ti)) 
is invariant under each of the companion matrices Hi(cri),..., An-d{cr\) (as a consequence of 
Lemma 1^ see equation (20)). This is perhaps the most crucial benefit of the regularization 
step that it makes the compatibility condition and the recursion independent of each other, 
plus it ensures that once the compatibility condition is satisfied, it is automatically satisfied in 
every subsequent step of the recursion. 


4 Discrete version of Noether’s normalization 

There are many autonomous n-D systems that are not strongly relevant. For example, the 
scalar system 53 = ker(criiT 2 — ui — cr 2 +1). Here the quotient module is not a finitely generated 
module over either M[(Ti, or M[cr 2 , cr^^]; neither is it a finite dimensional vector space over M. 
Therefore, this 53 is not a strongly relevant system of any order. In fact, autonomous systems 
are more often not strongly relevant. For such systems, clearly, the method of solution presented 
above becomes inapplicable. However, we shall see in this paper that by a normalization process 
every autonomous system can be rendered to an equivalent strongly relevant system. We call 
this normalization process the discrete Noether’s normalization lemma (DNNL); for it is an 
analogous version (for Laurent polynomial rings) to the well-known Noether’s normalization 
lemma (for polynomial rings). The result was discovered by E. Noether (1882-1935), who 
first presented it in her 1926 paper, entitled “Der Endlichkeitsatz der Invarianten endlicher 
linearer Gruppen der Charakteristik p”, and used it for studying invariant theory of finite 
groups over fields of arbitrary characteristics. For detailed analysis and applications of this 
result see |Eis95j . 
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Like the classical Noether’s normalization, our DNNL, too, crucially involves a coordinate 
transformation. However, the space on which we need to do the coordinate transformation 
differs because we are dealing with Laurent polynomial rings. It turns out that for DNNL we 
need to do a coordinate transformation on the indexing set Z”". 


4.1 Coordinate transformations in Z” and their effects on A and 


The idea of change of coordinates and its effects on a behavior are not new (see [Zam98|IVaini| ). 
however, we use this idea to a completely new end. By a coordinate change we mean a Z-linear 
map from Z"" to itself of the form 


T -.7/^ Z’^ 

C0l(l/i,Z/2, . . . jt'n) =: Tv, 


where T G ig a matrix whose determinant is ±1. Such matrices are called 7-unimodular. 
Note that det(r) = ±1 ensures that the columns of T span the whole of Z”" as a Z-module. 
Such a coordinate transformation T induces the following M-algebra homomorphism from A 
to itself. 


ifT '■ A 


.4 


(38) 


for all G Z"’. 

Since T is Z-unimodular, (pT turns out to be bijective] that is, pT is an automorphism of 
the M-algebra A. As a consequence, an ideal a C ^ is mapped to another ideal ptA)- The 
map pt can be extended to a map from A^ to itself: 


tpT • A^ —^ A^^ 

[/l(o-) /2(o-) ••• fqA)] [ttUiA)) TTU2A)) ••• ifTifqA))]. 


(39) 


The map pT is an .4,-module morphism via the automorphism i.e., for r{a) G A^ and 
/(cr) G A, 

^r(/(cr)r(cr)) = (^T(/(o-))^r(r(cr)). 

The bijective property of (pT extends to the module case: as a result, (Pt{TI), the image of a 
submodule TZ T A^ under (pT, too, is a submodule. 


4.2 Transforming a single polynomial 


We start off by normalizing a single polynomial. In Lemma 13 below we show that given 
any Laurent polynomial / there exists a coordinate transformation matrix T such that the 
image of / under the corresponding automorphism px has a special structure: pxif) is a 
linear combination of distinct powers of an with all coefficients being units in the smaller ring 
An-i := ..., cr„_i, This result is not new; a proof of this result is given in 

|Par04) ■ where it was used in the context of designing inverse filters for n-D discrete filters. 


We provide a proof of Lemma 13 here to make this exposition self-contained. 


Lemma 13. Let f € A be given by 




, where G M, and only finitely many 0. 




Then there exists a co-ordinate transformation T G such that under the corresponding 

homomorphism px '■ A ^ A, we have 

TtA) = ('^Uia^^ u, (40) 


,i=0 


where Ui G An-i and u G cr„^] are units and k is some finite positive integer. 
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Proof : In general, there are many T’s which will render / in the form of equation (40); we 
construct one particular T for which the monomials of / are mapped to monomials having 
different Un powers. To this end, let us define T G as 


1 0 ••• 0 

0 1 ••• 0 

ti t2 ■ ■ ■ 1 


where ■ ■ ■, tn-i G Clearly, T is unimodular. The homomorphism (px corresponding to 

this T maps a monomial a'' with u = col(z^i, 1/2,..., r'n) G Z"’ to the following 

= <1(72^2 ... 


Let us define t := col(fi,t 2 , • • ■ Tn-i) £ and t := col(t, 1 ) G Z"". Note that the power 

of an in ipT{cr^) is then given by the usual dot product (t, v). Let us also dehne the following 
finite subset of Z"'. 

supp(/) := {i^ \ 0}. 

We now claim that t can be chosen such that the cj,i-degrees of ipT{cr^) for v G supp(/) are all 
different from each other, i.e., the following set 


|t G Z^-^ 


(t, v) / (t, v') for all u, v' G supp(/), v 7^ 7^ 0 . 


(41) 


To verify this claim let us choose two distinct elements from supp(/), say u = col(z/i, 1 / 2 ,..., Vn) 
and u' = col(i/(, ..., and consider the following equation 


ui - 

V2 - V 2 


L Vn-\ - W-i J 


(42) 


We consider solutions to this equation (42) with t G 


»n—1 


I.e., 


Su,u' — {t G 


y>n—l 


equation (42) is satished} . 


The solution set if nonempty, is an affine hyperplane of dimension (n — 2). This follows 

from the following argument. Since u 7 ^ v', the vector col((z/i — z^(), ( 1^2 — 2 ^ 2)5 ■ • • > — 

and the number — Un) cannot both be equal to zero simultaneously. Hence, S^y is either 
empty or an affine hyperplane of dimension (n — 2). Now consider all pairs I'ji'' G supp(/) 
and the corresponding Si,y^s. Clearly, there are only finitely many such pairs (in fact, their 
number is ). Consider the following union 


5 := 


U 


z/'E supp (/), 1/7^ z/' 


Now observe that the set under consideration in equation (41) is related with the above defined 
S as 


t G Z' 


n—1 


(t, u) 7^ (t, u') for all u, u' G supp(/), z/ 7^ z/'| = (M” ^\ 5 ) PI Z' 


n-l 


Since the number of pairs {v, v') is finite and each of S^y is either empty or a hyperplane of 
dimension (n— 2 ), the union S cannot contain all of And therefore, (M”'“^\5)nZ"'“^ 7 ^ 0, 


which proves our claim equation (41) 
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It then follows that we can choose a t G such that 

(t, v) 7^ (t, v') for all u, v' G supp(/), v / u'. 

In that case, the the nonzero monomials appearing in (prif) will have all distinct powers of 
an- Therefore, ^rif) can be written as 



V>Tif) = Y1 

vesupp(/) 

Note that 0^=1 units in An-\- Rearranged in the ascending order of powers of 

an we get Uj’s from Wj=i Further, taking out the smallest power of an and calling 

it u we get the required expression of equation (40). □ 


There is an important consequence of Lemma we state this as Lemma 14 below. This 
will be essential for the general case Theorem |17| where we deal with normalization of a general 
(non-principal) ideal. 

Lemma 14. Let {0} 7 ^ o C ^ 6 e an ideal. Then there exists a coordinate transformation 
matrix T G inducing (pT '■ A ^ A, such that 

• A/^pt{o) contains An-i/ ^ An-i) as a subring, and 

• A/<pT{if) is a finitely generated module over An-i/ ^ An-i)■ 

Proof : First note that for an ideal a C the quotient ring Aja. can be written as 

A / An— I r 

■ 

a n An-I 

This is the algebra generated by On and over .A„_i/(an An-i). As we have already 

seen, ipT takes an ideal 0 C ^ to an ideal b = T A. Suppose / G 0 is any nonzero 

Laurent polynomial. By Lemma 13 there exists a coordinate transformation matrix T such 
that p>T{f) has the form given by equation (40). Since u in the right-hand-side of equation 
(40) is a unit in A, it follows that 


Trif) = J u G b 


G b. 


^i=0 


i=0 


Further, since {ui} are all units, the coefficient of the highest power of an, i-c. Uk, can be 
made 1, and the resulting polynomial will still belong to the same ideal. That is. 


Ufe (Tn 


k-l 

i=0 




G b 44 a. 


k-l 

k I \ ^ 


i=0 


where u[ = UiUj^^ G An-i- In other words, an satisfies a monic polynomial equation in A/h 
with coefficients coming from An-i/{fo H An-i)'- 


k-l 

+ ^u[{any = 0. (43) 

i=0 

Now, suppose we have a Laurent polynomial, say g £ A, which can be written as a poly¬ 
nomial with non-negative powers in an with coefficients from An-i- That is, g G An-i\an\. 
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Then, we can carry out Euclidean division algorithm by the monic polynomial appearing 
equation (43) above to obtain 


m 


9 = \(^r. 


k-l \ k-1 

i=0 / i=0 


where q £ An-ilcrn] and {ao, Oil, ■ ■ ■, «fc-i} ^ ^n-i- Equation (43) then implies 


fc-i 


9 


= '^ai {anY 
1=0 


Note that the coefficients above are all in An-i/{^ n^n-i), and g G ^ i \An]■ Thus, every 
element in bn.A~^i written as a finite linear combination of {l,a^,..., 

with coefficients from An-i/{h H An-i)- In other words, [oT] is an ^ i )-aiodule 

ra.t.prl bv -f 1 TfZ ( (T_ 


generated by {1, ..., ^}. 

Further, multiplying both sides of equation (43) by we get 


fe-i 

+ ^u'- {Ony~^ + u'o (^)“^ 

i=l 


= 0 . 


(44) 

i=l 

Since Uq was obtained by multiplying two units, namely uq and it follows that Uq is a unit 
itself. This is because units in A go to units in A/h. Therefore equation (44) can be rewritten 
as 

k-l 

= - (m'o) ^ - X] (45) 

i=l 

By taking positive powers in both sides of equation j4^, we get that for every f G N, the 
monomial (oT)~* can be written as a polynomial in bnX^t^]' Therefore, every polynomial 
in bn^~^ 1 [^~^] written as a polynomial in brf^~\ I^^ other words. 


An—l r_ It _ An—1 r_■ 


b n An-i 


[CTn 


b n An—i 


an]. 


Hence 


An—l r __It An—1 r _T 

[an, an \ — . < [an\. 


(46) 


A/b = - -^- 

b n An—l ^ j_ 

We have already shown that [^T] is a finitely generated ^ j,a4~^ i ) ~™odule. Therefore, 

by equation (^46) above, we get that Ajb too is a finitely generated ^ bnA~^ i )'Module. Also 
note that ^/bcontains .A„_i/(b n An-i) as a subring. This completes the proof. □ 


Remark 15. In the situation of Lemma 14, it follows that every element in A/ (Pt{(i) satisfies 
a monic polynomial equation over the subring An-i/ H An-i) (see |AM691 Proposition 
5.1]). In this case A/(Pt{<Y is said to be an integral extension over the subring An-i/ H 
An-i), or equivalently, the map An-i/ {^T{o)l~^An-i) A/ipT{o) is injective and integral. The 
converse of this, that is injective and integral implying the bigger ring being finitely generated 
module over the smaller ring, is not true in general. However, if the bigger ring is a finitely 
generated algebra over the smaller ring, then the converse is also true f |AM69l Corollary 5.2]). 
In our case, for every 1 ^ d ^ n and an arbitrary ideal a C .A, the quotient ring A,/a is 
the algebra generated by {a^+i, (ctTkT)”^ , • • • (^)~^} over Ad/{o. n Ad). Thus, A/a is a 

finitely generated algebra over Ad/{ar\ Ad). Therefore, in all the cases that follow from now, 
the following two statements are equivalent: 
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1 . A ring map Ad A/a is injective and integral. 

2 . The bigger ring A/a is a finitely generated faithful module over the smaller ring Ad- 
We use the two situations interchangeably in the sequel. 


4.3 The general theorem 


Lemma above shows that for a nonzero ideal a, there exists a coordinate transformation 
matrix T G such that An-i/ —)■ A/ipT{(^) is injective and integral. Note 

that the smaller ring An-i/ n.4.„_i) may not be a Laurent polynomial ring because the 
intersection ideal n . 4^-1 may not be the zero ideal. However, it is desirable to have 


the smaller ring to be a Laurent polynomial ring itself. Theorem 17 shows how this can be 


done with a suitable coordinate transformation T. We shall require the following Lemma 16 to 


prove Theorem 17 As done in |KR00j for the conventional Noether’s Normalization Lemma, 
we also follow an induction strategy for the proof of Theorem dzi While Lemma [T4| provides 
the first step of the induction, the inductive step is what Lemma deals with. The lemma 
tells us that if we already have an injective and integral ring map Ad/{a H Ad) —>• A/a such 
that a n Ad A {0}, then by using a coordinate transformation T we can go down one step 
further and have Ad-i /n Ad-i) —)• A/ipT{A) injective and integral. 

Lemma 16. Let {0} A ^ A A be an ideal and 1 A d A n an integer such that Ad/ (a H Ad) 
.4./0 is injective and integral. Suppose that aCi Ad A { 0 }- Then there exists T G , inducing 
(fT ■ A 3 cr’^ cr'^'^ G A for all v G HA, such that Ad-i/{^PT{A)r\Ad-i) -3- A/ipt{A) is injective 
and integral. 


Proof : We first define a := a n Ad an ideal in Ad, and then apply Lemma |14| to a and Ad- 
Since a 7 ^ {0}, by Lemma 14, there exists a unimodular T G inducing the automorphism 
: Ad -3 Ad such that 

Ad-i/{^f(a) n Ad-i) -3 Ad/A>f(^) (47) 


is injective and integral. 

With this T we now define our desired coordinate transformation on Z” as T := 


T 0 

0 Iji-d 

HA^n_ Note that the restriction of T on the first d coordinate axes is T and that on the rest of 
the coordinate axes is identity. Therefore, the corresponding ipT : A A maps elements from 
Ad to itself. In fact, frlAd — Pf' Since ipj; is an automorphism of Ad, we get that 


PriAd) = Pf{Ad) = Ad- 


(48) 


Now, since action of ipx restricted to Ad matches that of ipj; on Ad we get <^^(a) = (^r(n)- 
Note that, since ipT is a bijection, we have ipT{(^ H Ad) = ^t{A) H ipT{Ad). So, we can write 

^Pt(a) = Pt{A) n ifriAd) = (pT{a) n Ad (by equation @). 

Since Ad-i A Ad, we get from the above equation 

(Pff/d) n Ad-i = <y7T(a) n Ad-i. 


Thus equation (47) can be rewritten as: 

Ad-i/{pT{A) n Ad-i) -3 Ad/{pT{Aj n .4^) 


(49) 


is injective and integral. 
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We now claim that Ad/{'-PT[<^) n Ad) —>■ A/^t{o) is injective and integral. In order to see 
this, first note that, like before, A/^pt{o) can be written as a finitely generated algebra over 
the smaller ring Ad/{^T{'^) H Ad) as 


A/ipT{o) = 


Ad 


(/?T(a) n Ad 


Wd+l, {cTd+i) 


-1 


) O’n, (<7n) ] 


So, in order to prove our claim, it is enough to show that each of the generators of the algebra 
written above is integral over the smaller ring. That is, for all d+1 ^ i ^ n,ai and {d^)~^ satisfy 
monic polynomial equations over Ad/i'^T{o.) Now, by assumption. A/a is integral over 

Ad /(a n Ad)- This means for d + 1 ^ z ^ n, and {a^)~^ satisfy monic polynomial equations 
over Ad/{a n Ad)- That is, for d + 1 ^ z ^ zz, there exist positive integers ki and r^, and 
Uio, an, ---, aiki-i G Ad and bio, hi,.--, hn-i G Ad such that we have 


ki I fc, — 1 I 

+ Uiki-ia-' + 


+ o-iiai + 0,0 


(u. ^ + • • • + 6,10. ^ + b, 


'iO 


G a 
G a 


(50) 


We then apply ipT to equation (50) and use the fact that for d+1 ^ z ^ zz we have (prio'i) = o, 
and (pTicr//^) = o“^. Thus we get 


+ '^Tiaiki-i)crt' ^ H-h ^pTiail)(Ji + ifriaio) 




(O^ + (/ZT(6ir,-l)(Oi ^ + • • • + <^T(6il)Oi ^ + ipTiho) 




-1 


G 

G 


ipxia) 

ipT{a) 


(51) 


As we have already seen, from the definition of (pT it follows that Ad is invariant under tpT, 
that is, ipT maps elements from Ad to itself. Therefore, tpxiaio), ipTian),..., ^riciiki-i) G Ad 
and (fTiho), </?t(6,i), . • •, ^pTihn-i) G Ad- Thus equation © shows that for all d +1 ^ z ^ zz, 
both Oi, (d7)~^ G A/<pT{h) satisfy monic polynomial equations with coefficients coming from 
Ad/i'PTio.) n Ad), which means A/^pT{h) is integral over Ad/ip’Tih) H Ad). That 

Ad/(<^T(a) n Ad) —^ A/ <pT{h) 

is injective is obvious. Therefore our claim that Ad/(<^T(ci) n Ad) —)• A/(/?T(fl) is injective and 
integral is true. 

By combining equation (49) with this we get the following two-stage tower of integral ring 
extensions: 

Ad-i/{<PT{h) n Ad-i) —> Ad/{^T{h) n Ad) —> A/<pt{o). 

It is a well-known result in commutative algebra that a finite tower of integral ring extensions 
is itself integral (see |AM69| Corollary 5.4]). Therefore, we get that 


Ad-i/(<^T(a) n Ad-i) —> A/ipt{o) 

is integral. Injectivity of the above map is clear. This completes the proof. 


□ 


We are now in a position to state and prove the main result: Theorem 17 We call Theorem 
[^the discrete version of Noether’s Normalization Lemma, or DNNL in short. 

Theorem 17. (DNNL) Let a C A 6e azz ideal. Then there exists a coordinate transformation 
T G inducing the ring map tpT ■ A 3 cr'^ cr'^'^ G A for all v G TA, and a positive 

integer 0 ^ d ^ zz such that the canonical ring map 

Ad A/>pT{h) 

is injective and integral. Equivalently, the following two conditions hold: 

• A/iprih) contains Ad as a subring, and 
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A/^pT{o) is a finitely generated module over Ad- 


Proof : As mentioned before, we give a proof by induction on the number [n — d). The first 
step of the proof follows from Lemma [14} We may assume that o {0}, because otherwise the 


theorem trivially holds with d = n. Since a {0}, by Lemma 14 it follows that there exists 
Ti : unimodular such that 

n A„_l) —^ A/ipTi{o) 


is injective and integral. The inductive step is given by Lemma 16 If ipTfia)r\An-i = {0} then 
the result of the theorem follows with d = n — 1. On the other hand, if O An-i A { 0 } 

then by Lemma 16 there exists unimodular T 2 G such that 


A1-2/{PT2{.PTfia)) n An-2) -t A/^PT2{PTfio.)) 


(52) 


is injective and integral. Note that ipT 2 ° PTi = PT 2 Ti- Thus equation (52) becomes 

An- 2 l{PT2T1 (a) n An-2) —t A/ ‘fT 2 Ti (a) 

is injective and integral. Next we look at ipT 2 Ti{o) n An- 2 - If this is the zero ideal then the 
theorem holds with d = n — 2. Otherwise, we once again invoke Lemma [TBl by which we have 
T 3 G Z”^"" unimodular, such that 

An-3/ {‘PT 3 T 2 T 1 (a) n An-s) -t A/ ‘-PT 3 T 2 T 1 (fl) 

is injective and integral. Then we look at ^T 3 T 2 Ti{o) O An-s and depending on whether this 
ideal is the zero ideal or not we invoke Lemma [T6| again or stop. We illustrate this process of 
recursively applying Lemma 16 in the following flow-chart Figure 


In the flow-chart, after the first step, if ^pt{A) O An-i A {0} then Lemma 16 is applied 
recursively, and in the process i is incremented at each step. This recursive process continues 
as long as <-Pt{A) O An-i A {0}- Note that the recursive process is guaranteed to stop after 
finitely many n) iterations because i is upper-bounded by n. When the process stops, we 
have the following situation 


An-i/{g:>T{A) n An-i) -t A/i-pt{o) ( 53 ) 

is injective and integral, plus v^t(ci) H An-i = {0}. This means we can rewrite equation ( [^ as 

A„-i A/(/?T(a) 

is injective and integral. Defining d := n — i we get our desired result. □ 


Remark 18. The number d is of special importance to us. This number is equal to the Krull 
dimension of the Laurent polynomial ring Ad- It is a known result in commutative algebra 
that the Krull dimensions of two commutative rings Ci and C 2 are equal if Ci C 2 is integral 


f |Eis95j L It then follows from Theorem 17 that the Krull dimension of A/gPT{A) too is equal 
to d. Yet another fact from commutative algebra says that Krull dimension remains invariant 
under isomorphism. Since ipT is an automorphism of A, A/a and A/g^T{A) turn out to be 
isomorphic to each other as M-algebras. Therefore, the Krull dimension of A/a is d. 


We illustrate the result of Theorem in the following simple example. 


-1 


1 -,^2-1^2 




Example 19. Consider the ideal a = {(Ji(T 2 — cti — <J 2 + ^) A M[(Ti, cr 
ring (J 2 , cr^^]/a is not integral over either M[(Ti, crj"^] or ]R[cj 2 , 

us that with a coordinate change the quotient ring can be made integral over a smaller Laurent 


The quotient 
tells 


Theorem 


17 
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polynomial ring. Note that dim(M[(Ti, crj” ,ct 2 ,ct^ ]/o) = 1. So the first reduction in the flow¬ 
chart of Figure yields the desired normalization. By Lemma 13 we can take T G to 

be 


T = 


1 0 
2 1 


The corresponding cpT turns out to be (px ■ cri o-i(T|,cr 2 e-)- < 72 . This px maps the generator 
of a to 


Pt{o'10-2 - (Ti - (72 + 1) = (TiCjI - CJlfjl - (72 + 1 = Ul ((72 - Crf - CT^ ^<72 -|- (7^ ^) . 

It then follows that px((i) = (o'i ~ ^2 ~ + crr^)- Note that the generator of px{o) 

is a monic polynomial in (72 with coefficients coming from M[(7i, (7j~^]. In other words, 
in M[(Ti, (T^^, (72, (T^^]/(^T(fl) satisfies a monic polynomial equation over M[(7i, (7 ^^]/((/?t(ci) H 
M[(Ti, (7^^]). Moreover, ipx{o) H M[(7i, crj”^] is clearly the zero ideal. Thus we get M[cji,(7^^] —)■ 
M[(7i, (7j"^, (72, (T^^]/(/5r(a) is injective and integral. 
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4.4 DNNL extended to modules 


In this subsection we show how the normalization process can be extended to modules. This 
extension is essential in order to apply the normalization process to non-scalar systems. The 
extension of DNNL to modules is obtained by two crucial observations, which we elaborate 
now. In order to proceed further, we require the following algebraic object: the annihilator 
ideal of a quotient module A^, which is defined as follows: 


ann(A^) = {/(<t) G A \ f{cr)m = 0 G for all m G A^}. (54) 

Note that A4 is naturall'^a module over A/ann(A4). Further, since the free module A'^ is of 
finite rank (= q) over A, it is hnitely generated as a module over A. Now, A4, being a quotient 
of A‘^, is clearly finitely generated over A. Hence, A4 is finitely generated over A./ann(A4) too. 
This is the first crucial observation. 

The second observation involves the maps induced by a coordinate transformation of Z”’. 
We write this observation as a proposition below. 


Proposition 20. Let T G be a Z-unimodular matrix. Further, let (pT and (pT he the two 

maps induced by T as in equations (38) and (39), respectively. Suppose TZ F A^ is a submodule 
of the free module A^. Then the following holds: 


(pT {ann{A‘^/TZ)) = ann{A^/(P t{TZ)). (55) 

Proof : Note that /(cr) G arm{A‘^/(P t{TZ)) if and only if for all r{cr) G A'^ we have f{a)r{a) G 
(fTiJZ). However, since (px is a bijection, for every r{a) G A'^ there exists a unique ri(cr) G A'^ 
such that (pT{ri{a)) = r(<T). Therefore, we can write /(<t) G am\{A'^/( pt{TZ)) if and only if 
for all ri{cr) G we have /(cr)^'r(n(o')) £ Tt{TZ). Recall that ipx is an automorphism of A. 
Therefore, there is a unique gicr) := (pf^{f{cr)). Then we can write 

/(cr)^'r(ri(cr)) = ipT{g{cT))tpT{ri{cT)) = ipTig{cT)ri{a-)) G ipriTZ) ^(cr) G ann(.) 

But g{(T) G ami{A^/TZ) if and only if /(cr) G (/?T(ann(A.'^/7?,)). □ 


We now present the non-scalar or module theoretic version of DNNL. 

Theorem 21. (Module theoretic DNNL) Let TZ F A'^ be a submodule. Suppose the 
Krull dimension of aiin(A'^/TZ) is d. Then there exists a coordinate transformation T G Z"'^"' 
inducing the M.-algehra homomorphism px ■ A3 c'' G A for all v G TA, and the corre¬ 

sponding A-module homomorphism (px '■ A'^ -3 such that the quotient module A^^/(pxiJZ) 
is a finitely generated module over Ad- 


Proof : Let T G Z"'^"' be an arbitrary Z-unimodular integer matrix. We have already seen that 
for any submodule 7Z oi A^, the quotient module A^ jLZ is finitely generated as a module over 
A/am\{A^/7Z). In particular, since pxiJZ) too is a submodule of we get that A^/px{TZ) is 
a finitely generated module over A/ann(A'^/^T(’^))- Now, let us define o := ann(A'^/7^). By 
Proposition 20, px{a) = ann(A'^/^T(’^))- Therefore, A^ jpx^jV) is in fact a finitely generated 
module over Ajpxia). 

Since Krull dimension of A/a is equal to d, by applying Theorem 17 to a we get that 
there exists a Z-unimodular T G Z”^*^ such that A/px{a) is a finitely generated module over 
Ad- Thus, by combining this with the fact that A^^/(pxiJZ) is a finitely generated module over 
A/px{a) we get the following chain of modules 


AVpxiTZ) 

t 

A/px{a) 

t 

Ad 


^‘Naturally’ here means the scalar multiplication for the base rings A and Gl/ann(AI) are identical. 
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where each step is finitely generated. Therefore, /(pt{T^) must be finitely generated over 
Ad- Indeed, suppose {gi,... ,gk} ^ A^^/(pTiJ^) is a generating set for A ^as a module 
over A/ipT{o)-, and {hi,..., h^} C A/ipT{o) is a generating set for as a module over 

Ad- Then 

\^gihj .£{ 

is a generating set for A ^over Ad- □ 


We illustrate this result in the following example. 


Example 22. Consider the submodule TZ of <72, cr 2 generated by the rows of the 

matrix 


R(a) 


ai-l 2 
1 0 - 2-1 


It can be checked that the the quotient module A 4 = M[o-i,o-j”^,0-2,o-^^]^/ 7 ^ is not finitely 
generated over M[o-i,o-^^] or M[o-i, o-j”^], neither it is a finite dimensional vector space over M. 
The annihilator ideal of the quotient module in this case turns out to be a = (o-io-2 —0-2—o-i — 1). 
Take the unimodular matrix T = [ 2 ? ] for coordinate transformation. Then we get that under 
the ring map ipT the ideal 0 is mapped to the ideal = (o-io-f — 0-2 — o-io-^ — 1). Note that 

M[o-i, o-^^, 0-2, o-^^]/(/ 5 r(a) is a finitely generated module over ]R[o-i, o-jf^]. 

The extension of (fx to M[o-i, 0 - 2 , o-^^]^, that is, (px maps TZ to the row-span of the 
following matrix: 


R{a) 


(Ticrl — 1 2 

1 0-2-1 


Note that the annihilator of M[cri, <72, <72 ^]‘^/px{TZ) is given by It then follows that 

M[fTi, CTi^, (72, a 2 ^]‘^/pxiR-) is finitely generated as a module over M[cJi, erf ^]. Indeed, it can be 


checked that, 



is a generating set for M[(7i,(7]^ ^,( 72 , 0-2 ^]‘^/px{TZ) as a module over M[(7i,(7^ ^]. 


5 Representation formula for general autonomous systems 


In this section we make use of Theorem 21 to obtain representation formula for general au¬ 
tonomous systems. This clearly is indicative of the fact that this representation formula that 
we are going to propose will involve a coordinate transformation on the indexing set Z”. There¬ 
fore, we first describe the effects of a coordinate transformation of Z"- on the behavior of an 
n-D system. 

Recall that a coordinate trasnformation of Z"- is given by a Z-unimodular matrix T E Z"-^” 
as 

T : Z^ -> Z" 
col(i/i, 1 ^ 2 ,.. •, t'n) =: Tu. 


Now note that this T induces the following M-linear map from 

l"\1 . 


to itself: 


: 


y ^ 

w{u) i-> w{Tv), 


(56) 


This map ‘hr is often called the pull-back of the coordinate transformation T. The pull-back 
‘hr is a bijection because T is Z-unimodular. In fact, it is easy to see that 


(h„^ = ‘hr-i. 
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Recall the M-algebra homomorphism ipT : A induced by T, and its extension to mod¬ 

ules (fT '■ A'^ X A^. The following lemma shows how (px and the pull-back are interrelated. 

Lemma 23. Let v,w € be related by v = Then for r{a) G A'^ we have 

r{a)v = ^T{}PT{'f'{(^))w). 

Proof : It is enough to prove for the scalar case: for v,w ^ and /(<t) g ^ we must have 

f(a)v = ^T{pT{f{o-))w). (57) 


This is because tpT is defined as pT applied element-wise, and is linear. Further, note that 
every Laurent polynomial is a finite M-linear combination of monomials. Therefore, once again 
by linearity of ‘hr, it is enough to prove equation (57) for A 9 f{<^) = where G Z"" is 
arbitrary. 

Now, v,w £ are assumed to be related by u = ^t{w)- First note that for v, u' G TP, 
the action of the monomial on v is given by the following 


= v{v' + v) = = w{T{A + v)) 




Since A was arbitrary we get for all z/ G Z” 


cr’^v = ^T{pT{cr’^)w)- 


This is what we wanted to show. 


(58) 

□ 


With this we now bring out the effect of T on the behavior 18 of an n-D system. We show 
that the image of 18 under the pull-back is also a behavior. Not only that, in fact, the 


equation modules of 18 and <h7’(18) are related by pT in a way as shown in Proposition 24 
below. 


Proposition 24. Let TZ C A^ be a submodule with behavior and let T G be 

T-unimodular. Then we have 


mn)) = ^{pT-An)). 


(59) 


Proof : Let w G be an arbitrary n-D trajectory and let r{a) G A'^ be an arbitrary 

row-vector of polynomial shift operators. Then, by Lemma 23, we have the following identity 

r{a) {^t{w)) = {pT{'i"{cr))w) . (60) 


Since <pT and both are bijective, it follows from equation (60) that w G 18(^t( 77.)) if and 
only if ^T^w) £ IB (77). In other words, 

(®(^)) = ®(^t(77)). 

Since T is invertible we can start the argument with T substituted by T~^, and thus we obtain 
the desired equation (59). □ 


Proposition 24 allows us to give the following definition. 

Definition 25. (T-isomorphism of behaviors) Suppose * 81,182 G . Then I 81 is said to 
be T-isomorphic to I 82 if there exists a T-unimodular matrix T G Z”^*^ such that 

*Bi = $t(*B2). 
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Proposition 24 then translates in the language of Definition 25 to the following: 

Corollary 26. Let G have their equation modules TZi,TZ 2 , respectively. Then *Bi 

is T-isomorphic to 532 if o,nd only if TZi = (prp-i{TZ 2 ), which in turn is true if and only if 
7^2 = 

Combining Corollary (or equivalently, Proposition [2^ with the module theoretic DNNL 
(Theorem |21[) we get the following important result. 

Theorem 27. Let ^ ^ he an autonomous behavior with equation module TZ C . Further, 
let the Krull dimension o/ann (.A'^/7^) be equal to d. Then IB is T-isomorphic to a strongly 
relevant autonomous system of order d. 


Proof: First note that, by Theorem 21, there exists a Z-unimodular T such that under 

the induced module map (pT '■ ^ we have A^^/(pt{TZ) a finitely generated module over 

Ad- Therefore, is a strongly relevant autonomous behavior of order d, by Definition 

By Corollary 26 53 is T-isomorphic to iB(^'r(^))- Thus IB is T-isomorphic to a strongly 
relevant autonomous system of order d. □ 


Theorem 27 enables us to obtain the following representations for a general autonomous 
system. These representations follow from the representations of strongly relevant systems 
presented in Theorem and Corollary 10 Recall the first order system of equations given 
in Table 1. The representations given below involve these equations. The next result is a 
consequence of Theorem 

Corollary 28. Suppose ^ ST is an autonomous system with equation module TZ. Let d he 
the Krull dimension of aim [A^/TZ). Then, there exist a positive integer r and the following 
matrices: 


■) An-d{cri) G Af^ 


. X(^i) G AT, 

• Ti(cri),^2(o-i) 

. c(^i) G at , 

such that 53 admits the following representation: 


53 = <^ u G 


(61) 


3 X G , w G such that equations in Table 1 

are satisfied, and v = ^t{w!) 

The next result is a consequence of Corollary pR] 

Corollary 29. Let ^ G he an autonomous system whose equation module is TZ F A'^ 
Suppose the Krull dimension of ami{A'^/TZ) is d. Then there exist 

• a positive integer r, 

• a matrix X (cri ) G AT, 

• a set of n — d square matrices Afia-fi),... ,An-d{o-i) G A^^'^, 

• a matrix C{cr\) G A^^'^ 

such that 53 admits the following latent variable representation: 


53 = 


V G 


3 xG 


and w G satisfying 


_ xTi)____ 

^d+lTr Ail^CTi) 
<^d+2lr — A2{(T\) 

O'nTr An—di(^l) 

~-CiaT 


0 

d 

0 


= 0, u = 


(62) 
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With these representations for general autonomous systems, we are now in a position to 
provide an algorithm for finiding solutions for general autonomous systems. The algorithm is 


exactly the same as Algorithm 11 with one extra stage at the beginning, where we normalize 
the given system to a T-isomorphic system that is strongly relevant. 

Algorithm 30. (Solution formula for general autonomous systems) 


(Noether) Normalization 


Input: An autonomous system = ker(i?(<T)), where R{(t) G A*^'^. 


Computation: 

Define TZ = rowspan(ii(cr)) C 
Obtain a = ann(A'^/7^). 

Using Theorem 17 hnd out a Z-unimodular T G such that under 

the induced map ipT '■ Awe have Ad A/^t{o) injective and 
integral. 


Output: 

T. 

TZ := (pt{TZ). 


Regularization 


Input: TZ Q A'^. 

Computation: 

Obtain a generating set G := {gi ,..., Qr} C A^jTZ. 

Obtain the matrix of relations X{a\) G A!^'^. 

Obtain the companion matrices Ai(<ti), ... ,An-d{<^i) £ (ensure 

that these matrices are invertible in Ad)- 
Obtain the output matrix C(cri) G -4.^^’’. 


Output: The matrices X{cri), Ai(cri),..., A„d(<Ti), and C'(cri). 


Recursion 


Input: The set of initial data: Solve X{a\)x = 0 for x G 

Computation: 

for n = (z/i, 1 / 2 , • ■ •, i^n) £ set 

w(iy) = (C{ai) Ai{ai)’^‘i+ix'^ {ui, Vd). 

end 


Output: The sequence w := 


Renormalization 


Input: The sequence w := {w{v)'\v^in 

Computation: 

for V = (z/i, 1 / 2 , • ■ •, ^'n) £ set 
w{u) = w{Tu). 

end 


Output: The sequence w := {t(;(z/)}jygzr, 


6 The set of allowable initial conditions 


In Algorithms 


11 


and 
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the initial conditions x G 


are required to satisfy the equation 


X(cri)x = 0. This is in contrast with the 1-D situation, where one can always choose x in 
such a way that the set of initial conditions is free (see |Wil83) l. In this section we analyze the 
set of allowable initial conditions and address certain important issues, namely, freeness of the 
initial conditions, and non-autonomy of the state space TX. 
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6.1 Freeness of the set of allowable initial conditions 


We have already seen that the representation of an autonomous syetem given by Theorem 
can be viewed as an (n — d)-D flow equation arising out of an (n — (i)-D first order dynamics 
on an infinite dimensional space. We start by identifying this infinite dimensional space. 
Suppose 7^ C is a submodule such that A'^ jTZ is a finitely generated module over Ad- We 
fix a generating set {gi,... ,gr} for A?jH- as an ^d-module. Recall equation Q, where we 
identified A?ITZ with a quotient module of A^ via the map i/j : A^ —)• A^/TZ. Let X{a-i) G 
be a matrix of relations of { 51 ,..., gr}- We have defined the following set of vector valued d-D 
trajectories: 

JT := |x G I X{ai)x = o} , (63) 


as the state space. 

It is clear that ^ contains all the allowable initial conditions. Note that ZX is in fact a 
vector valued d-D behavior having X{(T\) as a kernel representation matrix. The submodule 
corresponding to this behavior is the row-span of the matrix X[a\) over Ad, which equals 
ker(^). Clearly, this submodule ]^er{^jJ) is unique for given set of generators of A'^/TZ as a 
module over Ad- Thus, by the one-to-one correspondence between behaviors and submodules, 
^ gets fixed once a generating set for A^/7Z (as an ^^-module) is fixed. However, even 
for a fixed TZ, different generating sets for A?jTZ as a ^^-module may exist. These different 
generating sets may even have different modules of relations, leading to distinct ZX?,. But, it 
can be shown that these different are isomorphic to each other in the sense of |PWR02] . 
because they have isomorphic quotient modules. One important question related to ZX is: 
when does A’^/7Z admit a generating set of cardinality r such that the corresponding ZX equals 
the whole of ? In the sequel, we refer to this situation by saying A^ jTZ admits a free 

X of rank r. 

It is easy to see that the answer to this question is: if and only if AZ^ jTZ is a free module of 
rank r over Ad- This is because then, and only then, A^ jTZ will admit a free generating set, 
i.e., the matrix of relations X(<ti) corresponding to the generating set will be the zero matrix. 
And consequently, X will be the whole of ) . This fact together with Swan’s extension 


|Swa78] of Serre conjecture to Laurent polynomial rings give us the following check (Theorem 


32) for determining whether AA jTZ admits a free X or not. 


Remark 31. In Theorem 32 we require the notion of a zero left-prime Laurent polynomial 
matrix. A Laurent polynomial matrix A(<ti) g A^^^ having full row rank is said to be 
zero left-prime if the ideal generated by the (k x k) minors of X(cti) is the whole ring Ad- 
See |Zer99l IBos82j for various different notions of prime-ness and their significances. One 


particular property of zero left-prime matrices that will be important for proving Theorem 32 
is the following: 

The row-span of a zero left-prime matrix is a direct-summand of the corresponding free module. 
And, conversely, if a submodule TZ of the free module A'^ is a direct-summand then TZ admits 
representation as the row-span of a matrix that is zero left-prime \Zer9tA \Bosh 


Theorem 32. Let TZ A'^ be a submodule such that A'^/TZ is a finitely generated module over 
Ad- Fix any generating set {gi,... ,gr} for A'^/TZ as an Ad-module. Then A^/TZ admits a 
free X of rank r' for some r' £ N if and only if the generators {gi,... ,gr} admit a matrix of 
relations X{a-i) G A//l^'' which is zero left-prime. 

Proof : Recall the map if : A/i —)• AA/TZ from equation (|^ and note that AA/TZ is isomorphic 
to ^0/ker(V') as ^d-modules. Now, as seen above, A^/TZ admits a free X of rank r' if and only 
if /TZ is a free ^^-module (of rank r'). By Swan’s extension of Serre’s conjecture to Laurent 
polynomial rings, A^/TZ is a free ^^-module if and only if A^/TZ is a projective ^^-module. 
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Since jTZ is isomorphic to ^JJ/ker(V') as ^d-modules, jTZ being projective is equivalent 
to ^Jj/ker(i/’) being projective. This, in turn, is equivalent to splitting of the following short 
exact sequence of ^^-modules ( lEisOSj i: 


0 —> keT{'ip) 


a:, 


•^d/ker(?/>) ^ 0. 


This is true if and only if ker(i/’) is a direct summand of A^^ ( |Eis95j ). But, a submodule of the 
free module A^^ is a direct summand if and only if it admits representation as the row-span 
of a zero left-prime matrix (see Remark 31). Thus, ker(^) being a direct summand of A^^ is 
equivalent to existence of an X(<ti) e A!^^ such that 


keic{'tp) = rowspan(X(<Ti)), 


where X[a\) is zero left-prime. 


□ 


Remark 33. The significance of having a free is that initial conditions can then be chosen 
arbitrarily. That is, solutions in ^(IZ) can then all be obtained by Algorithm 1 11 1 with arbitrary 


X E 


Note that in Theorem 
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it has been assumed that the behavior is strongly rele¬ 


vant. Therefore, in order to check freeness of initial conditions for a general system, one must 
first apply a suitable Noether’s normalization, and then invoke Theorem 32 This brings up 
the question whether freeness of initial condition is dependent on the particular normalization 
used, or is it an inherent property of the original system. This remains an interesting open 
question. 


6.2 Non-autonomy of the set of allowable initial conditions 

As mentioned above, the set of allowable initial conditions, or the state-space, ^ can be viewed 
as a d-D behavior itself. It is often desirable to have this d to be as small as possible. We show 
in this subsection that this question of whether d is the smallest possible or not is related with 
being non-autonomous as a d-D behavior. 

It is already well-known that an n-D behavior 53 is non-autonomous if and only if the 
quotient module A?ITZ is faithfu^ over A. Translated to the case of ^ it means ^ is non- 
autonomous if and only if A^/(pt{TZ) is faithful as an ^^-module. Keeping this observation in 
mind, it is then desirable to design the coordinate transformation T so as to ensure A'^/(pTiJZ) 
is faithful as an .A^-module. It so happens that our proof of the module theoretic DNNL 
(Theorem [ 2 ^ already guarantees this property for A^/(pt{TZ). We bring out this fact in the 
following Proposition which is a slightly modified version of Theorem \T\\ 

Proposition 34. LetTZ C A‘^ be a submodule. Suppose the Krull dimension of Qim{A^/7Z) is d. 
Then there exists a coordinate transformation T E inducing the "^.-algebra homomorphism 

ifT : A 3 (t‘' E A for all v ^ TA, and the corresponding A-module homomorphism 

(fT ■ A^ ^ A'^, such that the quotient module A^/ifTiTZ) is a finitely generated and faithful 
module over Ad- 


Proof : Recall the proof of Theorem 


We have already seen there that a unimodular 


T E that makes A^/(pt{TZ) a finitely generated module over Ad can be constructed 

using Theorem 17 a unimodular T E Z"^"" which renders .4./</?'r(ann(yl'?/77)) finitely gen¬ 
erated module over Ad suffices. However, the T constructed in Theorem does more; it 
makes the canonical map Ad —>• .4/<^'r(ann(^/77.)) injective and integral. This means the T 
constructed in Theorem 17 not only makes A/ipT{8x\^{A/7Z)) finitely generated over Ad, but 
faithful over Ad, too. In other words, by Theorem 17, there exists T E unimodular 


module M over A is said to be faithful if ami(AI) = {0}. 
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such that ^/y?T(ann(^/7^)) is finitely generated and faithful as a module over Ad- However, 
A'^/ipTi'J^) is clearly a faithful module over A/ipT{sinn{A'^ /TZ)). Thus A'^/ipTiTZ) is faithful 
over Ad- That A^^/pt^JZ) is finitely generated over Ad has already been proved in Theorem 

EH □ 


Combining this observation with the fact that ^ is non-autonomous if and only if A'^ /<pt(JZ) 
is faithful over Ad, we get the following interesting result. 

Corollary 35. Let ^ ^ be an autonomous n-D system such that the Krull dimension of 
A/aim{A‘^/7Z) is equal to d- Then there exists a unimodular T G such that 

1. is T-isomorphic to a strongly relevant autonomous system ^2 of order d, and 

2- the state space SL corresponding to 552 is non-autonomous as a d-D behavior. 

We have claimed in the beginning of this subsection that non-autonomy of is related 
with the smallest possible d. In order to prove this claim we first show that if the Krull 
dimension of .4./ann(^'^/7^) is equal to d then A'^/Zpt{'TZ) cannot be finitely generated as a 
module over Ad' for any d' < d- In other words, for the existence of a d'-T) behavior Z%Z it is 
necessary that d' ^ d. 

Proposition 36. Let TZ C A'^ be a submodule- Suppose the Krull dimension of A/axm{A‘^/LZ) 
is equal to d- Then there exists a unimodular T G A^/pt{TZ) is finitely 

generated as a module over Ad' if and only if d' ^ d- 

Proof : (If) We assume that d' ^ d, and we want to show that there exists a unimodular 
T G such that A'^/(pt{TZ) is hnitely generated as a module over Ad'- This is an imme¬ 


diate consequence of Theorem]^ (or Proposition 34). Indeed, by Theorem 21 there exists a 
unimodular T G Z”^”’ such that A^^/(Pt{TZ) is finitely generated as a module over Ad- However, 
since d' ^ d means Ad is a subring of Ad', A^^/priTZ) clearly is finitely generated as a module 
over Ad', too. 

(Only if) We assume that there exists a unimodular T G Z”^” such that A^jZpT^jV) is 
finitely generated as a module over Ad', and we want to show that d' ^ d- Let us define 


a := aim{A^/LZ)- Recall from Proposition 20 that (/^^(ci) = a.mi.{A‘^/( pt{TZ))- Since A'^/(pt{TZ) 


is finitely generated as a module over Ad', we get that A//(pT{TZ))a is also hnitely 


generated as an ^^z-module. Therefore, by the identity of Proposition 20, A/pt{(^) is hnitely 
generated as a module over Ad'- Therefore, the canonical map Ad'/ p>t{oZ) n Ad' A/pt{o) 


is injective integral (see Remark 15). Thus Krull dimension of A/ipt{<^) is equal to that of 
Ad' /pt{o) n Ad', but Krull dimension of Ad' /pT{of} H Ad' is clearly less than or equal to that 
of Ad', which is d'- Again, Krull dimension of A/ipT{(^) is d because A/pt{oZ) and A/a are 
isomorphic (see Remark 15). Therefore, d' d. □ 


Next we show that \i d' > d the Krull dimension of .4./ann(^‘i/7^), then the corresponding 
cannot be non-autonomous. 

Lemma 37. Let TZ C A^ be a submodule- Suppose the Krull dimension of Aj8xm[A^ jLZ) is 
equal to d- Further, let T G Z”^"’ be unimodular such that A^^/Zpt{TZ) is finitely generated as 
a module over Ad' - If d' > d then A^ cannot be faithful as a module over Ad' - 

Proof : Suppose, contrary to the claim, that A^/(pT^JZ) is faithful as a module over Ad'- Then 
the annihilator of A^jpT^JZfi as a module over Ad', which is ann(^'^/^'r(’^)) FAd', must be 


{0}. Utilizing the identity of Proposition 20 we get 


(^r(ann(^'^/7^)) Z^Ad' = {0}. 
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However, this means that the canonical map Ad' ^/y?'r(ann(^'?/7^)) injective. Further, 
as has been argued in the proof of Proposition 36, ^/</?T(ann(^'?/7^)) is finitely generated 
over Ad' because A^/(pTiTV) is so. Therefore, the canonical map Ad' —>• ^/<^T(ann(.A'?/7^)) is 
integral, too. Therefore, Krull dimension of A/(pT{3iim.[A^/TZ)) must be equal to that of Ad', 
which happens to be d!. This clearly contradics the assumption that d! > d. □ 


These observations can be summed up as follows: 


Theorem 38. Let ^ be an autonomous behavior with equation module TZ <£ A'^. Let d 

be the Krull dimension of A/aim{A‘^/7Z). Then there exists a unimodularT £ such that 

A'^/(pT^JZ) is finitely generated over Ad' with the corresponding state space XZ non-autonomous 
if and only if d' is equal to d. 


Proof : (If) This implication is just a restatement of Corollary ^ 


(Only if) By Proposition 36, d' ^ d; for in that case no unimodular T £ exists such 


that A^I(pt(JZ) is finitely generated over Ad'- On the other hand, if d' > d then, by Lemma 37 


A'^/(pTiJZ) cannot be faithful over Ad', which means ^ cannot be non-autonomous. Therefore, 
d' = d. □ 


7 Concluding remarks 

In this paper, we presented novel representation formulae for solutions of discrete n-D au¬ 
tonomous systems. Our analysis is based on a crucial algebraic property: the quotient module 
being finitely generated over a proper subring of the ring of shift operators, A. Quite ex¬ 
pectedly, many systems do not have this property. Therefore, we first consider those systems 
whose quotient modules are indeed finitely generated over a subring Ad- We called such sys¬ 
tems strongly relevant of order d. We then showed that such systems admit a state space 
like first order evolutionary equations. Using this we provided an explicit solution formula for 
trajectories in such systems. Next, we considered systems that are not strongly relevant. For 
these systems we showed that they can be converted to strongly relevant ones by a suitable 
change of coordinates on the indexing set The crucial fact that enables this conversion 
is what we have called the discrete version of Noether’s normalization lemma (DNNL). We 
showed in this paper that given any submodule TZ C A'^, there exists a coordinated trans¬ 
formation T ■. TP ^ TP such that, under the corresponding module homomorphism fixiTZ), 
the quotient module A^^/tpTfTZ) becomes finitely generated as a module over a subring Ad- 
Using this DNNL we have shown how a general discrete n-D autonomous system can be con¬ 
verted to a strongly relevant system. Then, making use of the earlier representation formula 
for solutions of strongly relevant systems, and the inverse of the coordinate transformation we 
obtained representation formula for general discrete n-D autonomous systems. 

An interesting feature of the representation formulae presented in this paper is that the 
solution trajectories are obtained by making (n—dj-D flow operators acting on initial conditions 
that are trajectories defined over T'^. This is very much akin to the state space solution formula 
of 1-D systems. However, there is one crucial deviation: unlike 1-D systems, here it may not 
always possible to chhose the initial conditions freely. In this paper, we provide necessary 
and sufficient conditions for these initial conditions to be free. Another interesting question 
concerning the set of allowable initial conditions is that whether it can be made smaller. This 
issue is related with the number d. We show in this paper that the smallest value that the 
number d can take is equal to the dimension of the characteristic variety of the system. 

Many important questions relevant to the results presented in this paper have been left 
unanswered. For example: can the representation formulae be extended to non-autonomous 
systems? Another relevant question is: what is the minimum number of state variables? Is 
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this number intrinsic to the system, or does it depend on the coordinate transformation? We 
have also not answered questions concerning how existing computer algebra packages can be 
used to compute the various matrices used in the representation formulae. 


Appendix 

Proof of Lemma |4] 


follows from elementary commutative algebra. However, for a direct proof of Lemma [39l the 
reader may refer to [PPT^ . where n = 2,d = 1 case of the result has been proved. The 
extension to the general n, d case is straightforward, and hence we skip it here. 

Lemma 39. The following are equivalent. 

1. ITZ is finitely generated over Ad- 

2. For every 1 ^ i ^ n — d there exists a polynomial Pi{cr) G ann(^'J/7^) of the form 

(^dU + + • • • + ai,i{(Ti)ad+i + ai,o{(Ti), (64) 

where Li is a finite positive integer, ajj((Ti) G Ad for alll ^ i ^ n — d and l^j^Li — 1, 
with aifl{cri) a unit for all 1 ^ i ^ n — d. 

3. For every \ ^ i ^ n — d there exists a polynomial pfic) ^ A of the form 

(^dli + ai,Li-i(cri)cr^|7^ H-h ai^i{ai)ad+i + ai,o(o-i), (65) 

where Li is a finite positive integer, aij(<Ti) G Ad for alll ^ i ^ n — d and l^j^Li — 1, 
with aifi{cri_) a unit for all 1 ^ i ^ n — d, such that 

Pi{o-)ej G TZ 


In order to prove Lemma we need the following algebraic result: Lemma 39 Lemma 39 


for all 1 ^ j ^q. 


With the help of Lemma we now prove Lemma 
Proof of Lemma^ Part (1) This assertion follows from repeated application of Lemma 
and the fact that the maps fij commute pairwise. 

Part (2) We first prove the result for the special case when q = 1. 

Case 1 (q = 1): In this case 7^ C ^ is an ideal, and hence, ann(^/7^) = TZ. Recall 
that Lemma deals with the situation when A/TZ\s finitely generated as a module over Ad. 
Therefore, by Lemma 39, for every 1 ^ i ^ n — d there exists a polynomial Pi{cr) G TZ having 
the following form 


(^d+i + o.i,F-ii^i)(^dfi ^ H-^ ai^i{cri)ad+i + ai,o(cri), 


where L, is a finite positive integer, aij{cri) G Ad for all 1 ^ ^ n — d and 1 ^ j ^ Li — 1, 

with aifi{ai) a unit in Ad for all 1 ^ z ^ n — d. 

Now consider a Laurent monomial where n € Z. Since Pi{cr) is monic in ad+i and its 
coefficient for the zeroth power of ad+i is a unit in Ad, we can carry out long division of a'f_^_- 
by pi{a) to get 


c^d+i = q{o-)pi{(T) + bi,L-i{(Ti)a^li ^ + bi^L-2{(Ti)a^li ^ + • • • + 6i,i(o-i)o-d+i + 6i,o(o-i) 


Li-2 


In other words, af_^- is congruent modulo TZ to an ^(^-linear combination of {1, + 

~ 1 4 


+ cr 


d-\-i 




35 








Notice that the above mentioned procedure can be carried out for all 1 ^ i ^ n —d. Hence it 
follows that a Laurent monomial of the form 0?=^! congruent modulo TZ to an ^^-linear 

combination of the following monomials: 

{ n—d 

X\^d+i 

i=i 

(The indices of the monomials in Q correspond to the parallelepiped, in the positive orthant 
of the integer grid with sides having lengths Li — 1, L 2 — 1,... ,Ln-d — 1 units. See 

Figure [^) Since every Laurent polynomial is an ^(^-linear combination of monomials of the 
form \YiZi follows now that every Laurent polynomial is congruent modulo TZ to an 

^rf-linear combination of monomials in Q. 


0-3 



0 ^ ^ T/i 1 


( 66 ) 


Figure 2: The parallelopied when n — d = 3 

Hence, the images of the elements in Q under the canonical map A -» A/TZ forms a 
generating set of as a module over Ad- We denote this set of images by G. We claim that 
under this generating set /Xj’s admit matrix representations that are all unimodular over Ad- 
We first show that has a matrix representation ^i(cri) that is unimodular. In order 
to do this we first put an ordering of the elements in Q. Consider the lexicographic ordering 
of monomials with < 7^+1 crd -\-2 ^ ^ <7n) and write down the elements of G in descending 

order. Let us name the elements of G, in this order, as {gi,g 2 , - - - ,gL}, where L = W^l 
Now, note that when the generators are ordered like this then they satisfy the following set of 
equations: 


Ati(ffi) = 52 


diigii) 

5l(5Li+l) 


—“i,o(<^i)5i “ ai,i(cri)52 — • • • — ai^Li-i{(^i)gLi 

9Li+2 


5i(52Li) 


-ai,o(<^i)5Li - ai,i(<^i)5Li+i - • • • - ai,Li-i(cri)52Li 


(67) 


It then follows that under this generating set (ordered in the above manner), admits a 
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matrix representation of the following form: 


Ai(cri) = 


^l.blockCcTl) 0 

0 ^l,block(o'l) 


where Ai,biock(o-i) G is given by 


0 

0 

^l,block(<^l). 


Gll,block(<7’l) = 


0 

0 


1 

0 


0 

1 


0 

0 






Note that det(yli^biock) = ai,o('^i)) a unit in Ad- Therefore, vli,block is unimodular over Ad- 
Hence, so is Ai{a-i). 

We can further claim that not only the map in fact, every one of the maps ^ 2 -, ^ 3 , • • •, fJ^n-d 
has matrix representations under the chosen generating set each of which are unimodular over 
Ad- In order to show this let us consider an i G {1, 2,..., n — d}- Now, we do the following 
trick: we reorder the elements of Q so that now once again they are in descending order with 
a different lexicographic ordering, namely, ad+i ■< (Jd+i -< (yd +2 ■< ■ ■ ■ < (Jn- Once again, we 
number the elements of Q in this order as {ffi, 52) • • • > d'h) where L = Li. Now, following 

exactly the same logic as above we obtain the fact that under this generating set (with this 
new ordering), the map ju-i admits a matrix representation of the form: 


Ai{cr-i_) = 


^i,block(*^l) 0 

b ^i,block(^l) 


where Ai_biock(<7-i) G is given by 


0 

0 


^i,block(<^l)J 


^i,block(*^l) — 


0 

0 


1 

0 


0 

1 


0 

0 






Therefore, because aj^o(<^i) is a unit in Ad, it follows that Glj biock(<^i) is unimodular over 
Ad- Hence, y4j(cri) is also unimodular over Ad- Now note that the new generating set 
{g'i,g 2 , - - -, 5^} is related with the old generating set {51,52, • • •, 5 l } by a permutation. Let Pi 
be the corresponding permutation matrix- That is, 


' 9 i 


’51’ 

92 

= P ^ 

52 

. 9 ' l . 


. 9 L _ 


Then, clearly, P~^Ai(ai)Pi =: is a matrix representation of the map gi under the old 

generating set { 51 , 52 , • • • ,5l}- The unimodularity of Ai{a-i) then easily follows from the fact 
that Ai[a\) is unimodular. 

Case 2 (general q): When 5 > 1 every element in A^^ can be written as 


ri(cr)ei + r2(cr)e2 H-h rq{o-)eq, 
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where ri(cr) G A for all i, and e^’s are the standard basis row-vectors of Recall that A^ITZ 
has been assumed to be a finitely generated module over Ad- By Statement 3 of Lemma 
then it follows that for every 1 ^ i < n — d there exists a polynomial Pi{cr) G ^ of the form 

CTd+i + H-h ai^i{ai)ad+i + ai,o(o-i), (68) 

where L, is a finite positive integer, ajj(<Ti) G Ad for all 1 ^ i ^ n — d and 1 ^ j ^ Lj — 1, 
with ai,o(f’'i) a unit for all 1 ^ z ^ n — d, such that 

Pi{(T)ej G TZ (69) 


for all 1 ^ j ^ q. Suppose we define, like Case 1, the set of monomials 

{ n—d 'I 


Then, by equations (68) and (69), it follows that for every 1 ^ i ^ q, ri{cr)ei is congruent 
modulo TZ to an ^^-linear combination of elements in Q times Cj. In other words, if we define 
the following set 

^vec := {mej | m G 1 ^ j ^ g} , 

then the images of under the canonical map A'^ -» A^ jTZ is a generating set for A^ jTZ as 
a module over Ad- Now, notice that, under this generating set, the map pi is represented by 
a matrix 


A ■ — 

'^2,vec — 


Ai 0 
0 Ai 


0 0 


0 

0 


where Ai is the representation of pi\ey By Case 1, Ai is unimodular. And hence, so is Ai 
Part (3) Note that if Ai((Ti)’s are invertible then the following diagram commutes: 


•^d 




Ai{a-i) i 


M 
7 Pi 
M 


-1 


(70) 


Repeated application of Lemma|^ along with equation (70) proves that the diagram of equation 
commutes. 

Part (4) This assertion follows from Part (3). 
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